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Global Impact Analyses 



9.1 INTRODUCTION 

The primary objective of this project was to determine the condition of the two World Trade Center 
(WTC) towers immediately following the aircraft impacts. This assessment includes both the 
determination of the structural damage that degraded their strength and the condition and position of 
nonstructural contents such as partitions, workstations, aircraft fuel, and other debris that influenced the 
behavior of the subsequent fires in the towers. The global impact analyses were the primary method by 
which most of the information on the tower damage was obtained. The analyses described in the other 
chapters of this report, such as the component and subassembly impact analyses, were performed to 
support the development and validation of these global impact models. 

9.2 ASSUMPTIONS AND LIMITATIONS 

A fundamental limitation of the global impact simulations was the maximum finite element model size 
that can be executed on the available 32-bit computer clusters. Since LS-DYNA uses addressable 
memory for domain decomposition, the combined aircraft and tower model could not exceed 
approximately 2.3 million nodes. The primary assumptions and limitations of the global impact analyses 
were the result of reducing the model size to meet these decomposition limitations as well as to achieve a 
run time that allowed the necessary global impact analyses to be completed within the duration of the 
investigation. 

Although the analyses were performed on a 32 bit computer cluster, the precision used in the analyses can 
be controlled by the analysis software. Both single precision and double precision versions of LS-DYNA 
were available for the impact analyses. In general, single precision analyses are more efficient and the 
precision is sufficient for the type of impact simulation being performed. However, when the dimensions 
of the structure being analyzed are sufficiently large, the single precision analyses can introduce rounding 
errors in the analyses. The rounding errors occur since the analysis is resolving deformations or 
analyzing element penetrations on a local scale that is several orders of magnitude smaller than the 
controlling dimension. 

In preliminary simulations, the coordinate system for the models of the tower structures was located near 
the base of the tower. As a result, the largest dimensions were the vertical position of the structures in the 
impact zone. This large vertical dimension controlled the size scale in the impact analyses and introduced 
rounding errors that were manifested as unstable element behaviors. To eliminate this precision problem, 
the tower model coordinate system was moved to be centered on the impact zone of the tower. The 
largest controlling dimension was therefore the distance across tower (significantly smaller than the 
height of the tower). After adjusting the coordinate location, the unstable element behaviors were no 
longer observed. 
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To confirm the adequacy of the single precision analysis, subassembly impact analyses (Chapter 6) were 
performed on the same model in both single and double precision. The comparison of the two analyses 
showed no substantial difference in the impact response and damage. 

Specific assumptions and limitations introduced in the analyses to meet the computational constraint 
included: 

• Tower structures away from the impact zone had a coarse mesh resolution. Therefore, 
damage in those regions may not be as accurately captured as that in the impact zone. An 
example is the potential damage to the exterior panels on the far side of the tower (opposite to 
impact). As debris passed through the building and struck a panel on the far side of the 
building, the coarse mesh refinement and merged boundary conditions at column ends (as 
opposed to bolted connections in the impact zone) underestimated the secondary impact 
damage. 

• Tower contents (workstations and nonstructural walls) were only included in the expected 
path of the aircraft impact and subsequent debris cloud. Therefore, debris and fuel that 
passed beyond this region could move more freely through the structure, only interacting with 
primary structural components. 

• The analysis of the impact response of the aircraft fuel cloud had several limitations. Smooth 
Particle Hydrodynamics (SPH) was used to model the fuel in the impacting aircraft, and 
approximately 60,000 SPH particles were used for the approximately 60,000 lb of fuel in 
each aircraft. Therefore, the particle size in the fuel cloud was approximately one pound. 
The air in and around the towers was not modeled, so the deceleration of the fuel particles in 
the cloud by aerodynamic resistance was not included. The contact algorithm for the fuel 
particles and tower did not include a sticking or "wetting" behavior so the fuel particles 
bounce off of components in the tower. The results of these limitations would spread the fuel 
cloud over a larger region in the simulation. Finally, the deflagration of the fuel was not 
modeled, and the resulting dynamic over-pressures in the tower from the combustion process 
were not included in the analysis. 

• Windows were not modeled on the exterior of the tower. The open space between the 
exterior columns allowed fuel particles and small debris fragments from the aircraft and 
tower to escape that may have been contained if the windows were included. Note that the 
weight of the windows was added to the columns as part of the superimposed dead loads 
(see Section 3.2.4). 

• The rotational velocity of the spinning aircraft engine components was not modeled. The 
effects of the rotational kinetic energy, spin stabilization of the engine trajectory, and 
potential for engine thrust during impact were, therefore, not included in the analysis. An 
analysis was performed to estimate the magnitude of this assumption, as described in 
Section 10.3. The analyses indicated that the rotational kinetic energy of each engine was 
approximately one percent of the aircraft initial kinetic energy. In addition, much of this 
rotational energy was probably dissipated by internal deformations of the engine components 
following impact with the tower exterior. Therefore, this approximation should have a small 
influence on the global impact damage. 
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• Aeroelastic forces were not applied to the aircraft wings since the generated stresses were not 
expected to affect the impact response. The energy associated with the aeroelastic stresses 
are much less than the energy absorbed during impact. Note that a wing tip deflection was 
applied as explained in Section 4.3. 

• Initial service loads (stresses) in the tower and aircraft were not included. The material 
internal energy associated with the elastic service loads were small compared to the material 
internal energy capacity. Therefore, it was not expected to have a significant influence on the 
dynamic impact response and deformation. Simplified analyses were performed to evaluate 
the magnitude of this assumption for the impact response of a core column, as described in 
Section 10.4. These analyses indicated that ignoring the static preload in the column had 
little influence on either the dynamic column deformations or the resulting collapse strength 
of the column. Gravitational acceleration was modeled during the impact analyses to include 
the gravitational acceleration on debris movement and potential contributions to the collapse 
of the damaged truss floor regions. 

• The impact analyses are subject to uncertainties in the input parameters, such as initial impact 
conditions, material properties and failure criteria, aircraft mass and stiffness properties, mass 
distribution inside the towers, the jet fuel distribution and dispersion, connections behavior, 
the presence of nonstructural building contents, etc. (see Chapter 8 for further details). The 
global analyses presented in this chapter not only provide a "base case" based on a best 
estimate of all input parameters (based on photographic and video evidence, material testing, 
and data in the open literature), but also provide more and less severe damage estimates based 
on variations of the most influential parameters, identified in the sensitivity analyses 
(Chapter 8). These more and less severe damage scenarios provide a range of damage 
estimates of the towers due to aircraft impact. 

9.3 MODELS FOR THE BASE CASE IMPACT ANALYSES 

The development of the tower models and the development of the Boeing 767 aircraft model were 
described in detail in Chapters 2 and 3, respectively. Similarly, the analysis of the aircraft impact 
conditions was described in Chapter 7. In this section, a brief description of the combined global model 
of the tower and aircraft used for the impact simulations is provided. Similar to the component and 
subassembly analyses, the global impact analyses were performed using the LS-DYNA general purpose 
finite element software package. 

The combined aircraft and tower models developed for the base case global impact conditions of WTC 1 
and WTC 2 are shown in Figure 9-1 and Figure 9-2, respectively. The global impact models were each 
in excess of two million nodes. A detailed summary of the number and types of elements in each model 
is provided in Table 9-1. The bulk of the aircraft and tower structures were modeled using shell 
elements. Solid elements were used for specific components, such as the exterior wall butt plates and for 
some of the thick hubs and flywheels in the aircraft engines. Beam elements were used primarily in the 
WTC tower floor truss structures and for the bolted connections in the impact zone of the exterior wall. 
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Table 9-1. Summary of the WTC global impact models. 




WTC 1 Tower Model 


WTC 2 Tower Model 


Number of Nodes 


2,068,736 


2,110,970 


Belytschko-Tsay Shell Elements 


1,682,615 


1,716,249 


Constant Stress Solid Elements 


73,189 


72,906 


Hughes-Liu Beam Elements 


47,952 


53,488 


SPH Fuel Particles 


60,672 


60,672 



The objective of the base case impact analyses was to provide an estimate of the impact scenario based on 
an analysis with a best estimate of all the variables obtained from photographic and video evidence, 
material testing, and data in the open literature. A summary of the aircraft impact conditions used for the 
base case global impacts of the towers is provided in Table 9-2. 



Table 9-2. Summary of base case impact conditions. 





WTC1 


WTC 2 


Impact Speed (mph) 


443 


546 


Vertical Approach Angle 
(Velocity vector) 


10.6° below horizontal 
(heading downward) 


6° below horizontal (heading 
downward) 


Lateral Approach Angle 
(Velocity vector) 


180.3° clockwise from 
Structure North 


13° clockwise from Structure 
North 


Vertical Fuselage Orientation 
Relative to Trajectory 


2° nose-up from vertical 
approach angle 


1° nose-up from vertical 
approach angle 


Lateral Fuselage Orientation 
Relative to Trajectory 


0° clockwise from lateral 
approach angle 


-3° clockwise from lateral 
approach angle 


Roll Angle (left wing downward) 


25° 


38° 


Total Aircraft Weight 


283,600 lb 


277,600 lb 


Aircraft Fuel Weight 


66,100 lb 


62,000 lb 


Aircraft Cargo Weight 


40,600 lb 


38,700 lb 
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(b) Side view 
Figure 9-1. WTC 1 global impact model. 
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(b) Side view 
Figure 9-2. WTC 2 global impact model. 
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WTC 1 BASE CASE IMPACT ANALYSIS 



This case is referred to as Case A for the remainder of the WTC Investigation reports. The WTC 1 base 
case analysis was performed for a 0.715 s duration following initial impact of the aircraft nose with the 
north exterior wall. The analysis was performed on a computer cluster using twelve 2.8 GHz Intel Xeon 
processors, each on a separate node of the cluster. The run time for this analysis was approximately two 
weeks. The progress of the global impact simulations was monitored on average every two days. The 
calculations were terminated when the damage to the towers reached a steady state and the motion of the 
debris was reduced to a level that was not expected to produce significant impact damage. A typical 
value for the residual kinetic energy in the airframe components at the termination of a global impact 
simulation was less than one percent of the initial kinetic energy at impact. 

A side view of the base case WTC 1 global impact model is shown in Figure 9-3. The tower interior 
contents were removed, and the tower structures were shown as transparent, so that the impact response in 
the tower interior could be seen. The floor numbers and core columns for the tower model are labeled in 
the figure. The corresponding side view was used to illustrate the base case global aircraft impact 
response in Figure 9-4. The response is shown at intervals of 0.1 s from impact through the initial 0.5 s 
of the response. Corresponding plan view images of the progression of aircraft impact response are 
shown in Figure 9-5. 
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Figure 9-3. WTC 1 base case global impact model. 

The initial 0.1 s of the base case global aircraft impact response, shown in Figure 9-4(b) and 
Figure 9-5(b), was dominated by the impact, penetration, and fragmentation of the forward fuselage 
structures. The engines and wing sections were just starting to impact the exterior wall. The forward 
fuselage structures were severely damaged both from the penetration through the exterior columns and the 
penetration of the 96th floor slab that sliced the fuselage structures in half. The downward trajectory of 
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the aircraft structures caused the airframe to collapse against the floor, and the subsequent debris motion 
was redirected inward along a more horizontal trajectory parallel to the floor. 
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(c) Time=0.20 s 
Figure 9-4. WTC 1 base case global impact analysis (side view). 
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(f) Time=0.50 s 
Figure 9-4. WTC 1 base case global impact analysis (side view) (continued). 
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(a) Time=0.00 s 
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(c) Time=0.20 s 
Figure 9-5. WTC 1 base case global impact analysis (plan view). 
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(d) Time=0.30 s 




(e) Time=0.40 s 
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(f) Time=0.50 s 
Figure 9-5. WTC 1 base case global impact analysis (plan view) (continued). 
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At 0.2 s after impact, the wings had completely penetrated the exterior wall, and only the tail structures 
were still outside of the tower, as shown in Figure 9-4(c) and Figure 9-5(c). The wing structures were 
completely fragmented by the penetration through the exterior wall. The aircraft fuel cloud was starting 
to spread out, but was still relatively dense, and the leading edge of the fuel was just reaching the tower 
core. By 0.2 s, the downward trajectory of the aircraft structures had transferred sufficient vertical load 
that the truss floor structures on the 95th and 96th floors were starting to collapse in the impact zone. 

At 0.3 s after impact, the aircraft was completely inside of the tower (full penetration completed at 
approximately 0.25 s), as shown in Figure 9-4(d) and Figure 9-5(d). The airframe was mostly broken up, 
but some large sections of the aft fuselage and tail were still intact, having penetrated through the opening 
in the tower wall produced by the forward fuselage structures. The aircraft fuel cloud penetrated 
approximately half the distance through the core and was spreading out. However, the subsequent motion 
of the aircraft fragments and fuel debris cloud began to be noticeably slowed beyond this time. The fuel 
and debris did continue to spread through the tower, but at a much slower rate, as seen in the remaining 
images in Figure 9-4 and Figure 9-5. 

The load transfer of the aircraft impact can be shown by the history of the resultant momentum of the 
aircraft as shown in Figure 9-6. The momentum plotted is for all of the aircraft structures and contents 
(including fuel) and is normalized by the initial momentum magnitude. The curve illustrates an initial 
rate of load transfer during the first 0.1 s of impact as the forward fuselage penetrated the exterior wall 
and impacted the interior structures. Between 0.1 s and 0.25 s, a more rapid load transfer rate was 
observed, as the area of the impact was greater (extending outward in the wing impact regions), and a 
higher percentage of the aircraft mass was impacting the interior structures. At 0.25 s, as indicated in 
Figure 9-6, the aircraft completely penetrated the building and retained approximately 30 percent of its 
initial momentum. Beyond this time, the rate of load transfer was steadily decreasing, with very little 
load transfer after approximately 0.5 s. 
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Figure 9-6. Normalized aircraft momentum for the WTC 1 base case impact. 
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The aircraft was severely broken into debris as a result of the impact with the tower, as shown in 
Figure 9-7. At the end of the impact analysis, the aircraft was broken into thousands of debris fragments 
of various size and mass, as shown in Figure 9-7(b). A closer inspection of the debris field shows that 
larger fragments still occurred for specific components such as the engines, shown in Figure 9-7(c). At 
the end of the simulation, the port engine was still inside the core, and the starboard engine was roughly 
one third of the distance from the core to the south exterior wall. Both had a speed of less than 50 mph. 
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(b) Aircraft debris field (time=0.715 s) 





(c) Major engine debris (time=0.715 s) 
Figure 9-7. Aircraft breakup and debris in the base case WTC 1 impact. 
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9.4.1 WTC 1 Base Case Impact Analysis - Structural Damage 

The structural damage to the WTC 1 tower by the base case impact conditions is described in this section. 
The primary structural components of interest were the exterior wall, core columns and core framing 
components, and the floor structures and concrete floor slab. The response of the internal contents, walls, 
and partitions were also of interest to provide the initial conditions for the fire dynamics simulations. 

Exterior Wall Damage 

The calculated damage to the exterior wall is significant for two reasons: (1) the exterior wall contributed 
significantly to the strength of the tower, and the degradation in exterior wall strength was important for 
the collapse analyses, and (2) the exterior wall was the one structural system for which direct visual 
evidence of the impact damage was available. Therefore, the comparison of the calculated and observed 
exterior wall damage can provide a partial validation of the analysis methodologies used in the global 
impact analyses. 

A comparison of the north exterior wall observed and calculated damage from the base case WTC 1 
global impact analysis is shown in Figure 9-8. The calculated impact damage to the exterior wall is 
shown with color fringes representing plastic strain magnitude, with undamaged sections in blue and 
strains at or above 5 percent shown in red. The schematic of observed damage was developed from 
inspections of the video and photographic data collected on the tower after impact. Both the observed and 
calculated damage wall regions shown in Figure 9-8 illustrate a region of the exterior wall from 
column 108 to column 152, extending from floor 91 to floor 100 (spandrels at floors 92 through 100). 

The first observation from the comparison of the calculated and observed damage was that the geometry 
and location of the impact damage zone were in good agreement. This agreement in the position and 
shape of the impact damage served to validate the geometry of the aircraft model, including the aircraft 
orientation, trajectory, and flight distortions of the wings. 

The second observation was the agreement in the magnitude and mode of impact damage on the exterior 
wall. The exterior wall completely failed in the regions of the fuselage, engine, and fuel-filled wing 
section impacts. Damage to the exterior wall was observed all the way out to the wing tips, but the 
exterior columns were not completely failed in the outer wing and vertical stabilizer impact regions. 
Failure of the exterior columns occurred both at the bolted connections between column ends and at 
various locations in the column, depending on the local severity of the impact load and the proximity of 
the bolted connection to the impact. The agreement of both the mode and magnitude of the impact 
damage served to partially validate the constitutive and damage modeling of the aircraft and exterior wall 
of the tower. 

Core Structural Damage 

The determination of the damage to the core columns and core beams was important in determining the 
residual strength for the subsequent analyses of structural stability and collapse. The overall model for 
the core structure and calculated impact damage to the core is shown in Figure 9-9. The figure shows 
that the core had significant damage in the region close to the impact point. The columns in line with the 
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aircraft fuselage failed on the impact side, and several of the core beams were also severely damaged or 
failed in the impact zone. 
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(a) Schematic of observed damage 
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(b) Calculated damage (t=0.715 s) 
Figure 9-8. Base case impact damage to the WTC 1 exterior wall. 



NISTNCSTAR 1-2B, WTC Investigation 



209 



Chapter 9 




(a) Calculated impact damage 
Figure 9-9. Base case impact damage to the WTC 1 core frame structures. 
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The calculated damage to the core columns by row is shown in Figure 9-10. The columns are shown 
with color fringes representing plastic strain magnitude with undamaged sections in blue and strains at or 
above 5 percent shown in red. A summary of the column damage is listed in Table 9-3 and shown 
graphically in Figure 9-11. The qualitative classification of the column damage levels is shown in 
Figure 9-12. The classification levels were light damage, moderate damage, heavy damage, and failed 
(severed). The light damage level was defined as having evidence of impact (low level plastic strains), 
but without significant structural deformations. The moderate damage level had visible local distortions 
of the column cross section (e.g. bending in a flange) but no lateral displacements of the column 
centerline. The heavy damage classification was for impacts that produced global deformation of the 
column, resulting in a permanent deflection of the column centerline. The severed columns were 
completely failed and could carry no residual load. A summary of the floor-by-floor damage to the core 
framing is shown in Figure 9-13. 

P-A effects generated due to the sway of the towers after impact, as observed in video evidence (see 
NIST NCSTAR 1-5A), were not expected to affect or impose additional damage to the core columns. 
The core columns were designed as axially loaded members without continuity of framing, and thus 
would not develop significant P-A moments (see Chapter 5 of NIST NCSTAR 1-2A). 

Floor Truss and Slab Damage 

The impact damage to the floor trusses and floor slabs could contribute significantly to the loss of 
structural strength and subsequent collapse. The truss floor sections provided lateral support to the 
exterior wall at each floor level. Any damage or holes in the concrete floor slab could provide a path for 
the fires to spread from floor-to-floor. Therefore, the calculation of the floor system damage was an 
important component of the global impact analyses. 

An overall front view of the floor truss structure in the impact zone, along with the calculated impact 
damage to the floor trusses, is shown in Figure 9-14. The truss structures shown are for the detailed truss 
floor section models described in Section 3.2.3 and Section 3.2.5. The figure is from a view point 
exterior to the tower, looking toward the impact point, and shows that the trusses experienced significant 
damage in the impact zone. A plan view of the calculated damage to the truss on each floor is shown in 
Figure 9-15. The calculated impact response produced severe damage to the truss structures in the 
primary impact path of the fuselage. The truss structures were severely damaged from the exterior wall to 
the core. The truss floor system on floors 94 through 96 were damaged and sagging downward as a result 
of the impact loading. 
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Figure 9-10. Base case impact damage to the WTC 1 core columns. 
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Table 9-3. Summary of core column damage for the base case WTC 1 impact. 


Column 


Location 


Damage Level 
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Figure 9-11. Summary of core column damage for the base case WTC 1 impact. 
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Figure 9-12. Classification of damage levels in core columns. 
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(c) Floor 95 core framing damage 
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(d) Floor 96 core framing damage 




(e) Floor 97 core framing damage (f) Floor 98 core framing damage 

Figure 9-13. Base case impact damage to the WTC 1 core floor framing (plan view). 
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(b) Calculated damage (t=0.715 s) 
Figure 9-14. Base case impact damage to the WTC 1 floor trusses (front view). 

The calculated damage to the WTC 1 floor slab for floors 94 through 97 is shown in Figure 9-16. The 
fringes of damage were set such that the concrete slab failed in the regions colored red (2 percent plastic 
strain was used corresponding to the zero strength strain limit for the concrete in unconfined 
compression). At these strain levels, the concrete was severely damaged and probably removed, exposing 
the supporting metal decking. Beyond 2 percent plastic strain, the strength of the floor slab was severely 
reduced in the analyses to model the residual strength of the metal deck after the concrete failure, 
breakup, and removal. At a plastic strain of 30 percent, corresponding to failure levels for the metal 
decking material, the elements were eroded (seen as holes ruptured in the floor slabs shown). The line of 
damage, visible around the edge of the core in Figure 9-16, is believed to be the result of numerical 
precision errors and interference in the modeling methodology at the connection of the floor slab rather 
than impact damage. 
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(b) Floor 95 truss damage 
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Figure 9-15. Base case impact damage to the WTC 1 floor trusses (plan view). 
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(c) Floor 96 slab damage (d) Floor 97 slab damage 

Figure 9-16. Base case impact damage to the WTC 1 floor slab (plan view). 

WTC 1 Base Case Impact Analysis - Fuel and Debris Distributions 



The objectives of the global impact analyses included the determination of both the structural damage and 
the initial conditions that influenced the initiation and propagation of the fires in the towers, including the 
distribution of the jet fuel in the towers, the distribution of tower contents and aircraft debris that provided 
flammable materials for the fires, and the condition of the partitions and walls that provided barriers to air 
flow and spreading of the fires. 

The global impact results presented in this section included the distribution of the jet fuel inside the tower, 
the damage to the building contents (partition walls and workstations), and the aircraft debris distribution 
in the towers. The jet fuel was modeled using the smoothed particle hydrodynamics (SPH) algorithm in 
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LS-DYNA, as described in Section 5.5.3. For the base case WTC 1 global impact analysis, the calculated 
distribution of the fuel in the tower and shape of the fuel cloud in a plan view and side view are shown in 
Figure 9-17 and Figure 9-18, respectively. At the termination of the global impact analysis, the residual 
momentum of the jet fuel in the impact direction was less than one percent of the initial momentum, 
indicating that the fuel cloud was nearly at rest at about 0.715 s. 

To more clearly present the calculated response of the structures that influenced the fire propagation, the 
structural components were removed from the visualization, with the exception that the core columns 
were maintained in the visualizations for reference positions. A plan view of the response of the 
remaining building contents and aircraft debris is shown in Figure 9-19. Similar plan views of floor-by- 
floor slices through the building contents and debris field are provided in Figure 9-20 through 
Figure 9-24 for floors 94 through 98, respectively. The bulk of the aircraft debris and fuel were arrested 
prior to exiting the far side of the tower core. A small amount of aircraft debris was calculated to exit the 
far side of the tower (south wall of WTC 1). 

A quantitative characterization of the fuel and aircraft debris distribution was obtained by slicing the 
model at vertical floor locations and calculating the mass at each floor level. A summary of the floor-by- 
floor fuel and debris distributions is given in Table 9-4. The bulk of the fuel and aircraft debris was 
deposited in floors 93 through 97, with the greatest concentration on floor 94. Approximately 18,000 lb, 
or 7 percent, of aircraft mass was eliminated from the debris cloud at the final state as a result of the 
erosion in the aircraft structures due to impact and breakup. This mass was not accounted for in the fuel 
and debris distributions provided in Table 9-4. A first approximation would be to increase the airframe 
debris distribution proportionately to account for the eroded mass. This eroded mass was maintained in 
the calculation but was no longer included in the contact algorithm. As a result, any residual momentum 
at the time of erosion could not be subsequently transferred to the tower. 

The calculated debris distribution included 17,400 lb of debris and 6,700 lb of aircraft fuel outside of the 
tower at the end of the impact analysis, either rebounding from the impact face (north wall) or passing 
through the tower (south wall). These estimates of mass outside the tower were expected to be 
overestimated in the calculation since the exterior walls were not modeled with windows that could 
contain the fuel cloud and small debris inside the towers. In addition, the impact behavior of the aircraft 
fuel cloud did not include the ability to stick to, or wet, interior components. Rather, the aircraft fuel SPH 
particles tended to bounce off of internal structures. 

The physics of fuel impact and dispersion in this type of impact event is complex and no appropriate 
validation data could be found. The fuel starts as a continuous fluid within the tanks and ends up 
distributed both on the tower structures and as small droplets that interact with the atmosphere 
surrounding the impact zone. No single analysis technique is currently available that can analyze this full 
range of fuel dispersion without significant uncertainties. 

Both the SPH and ALE analysis techniques (see Section 5.5.3) available for the analysis of the fuel 
impact and dispersion had limitations. Details of the fuel behavior such as the wetting of the fuel against 
tower structures and interior contents or the physics of the fuel breakup into droplets are not accurately 
reproduced in either analysis technique. However, the momentum transfer from the fuel to the tower 
structures and subsequent impact damage produced by the fuel can be modeled by both analysis 
techniques. 
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(a) Calculated impact response (floor slab removed) 
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(b) SPH fuel distribution 
Figure 9-17. Calculated fuel distribution in the base case WTC 1 analysis (plan view). 
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(a) Calculated impact response 




(b) SPH fuel distribution 
Figure 9-18. Calculated fuel distribution in the base case WTC 1 analysis (side view). 
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(a) Impact configuration (plan view) 
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(b) Calculated impact response (t=0.715) 
Figure 9-19. Calculated WTC 1 building contents, fuel, and aircraft debris distribution. 
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(c) Calculated impact response (fuel removed) 




(d) Calculated aircraft debris 
Figure 9-20. Calculated floor 94 contents, fuel, and aircraft debris distribution. 
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(a) Initial impact configuration 
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(b) Calculated impact response 







mn 

JUTU 



j li ii n 
(c) Calculated impact response (fuel removed) 



TPP-5 






(d) Calculated aircraft debris 
Figure 9-21. Calculated floor 95 contents, fuel, and aircraft debris distribution. 
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(b) Calculated impact response 
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(c) Calculated impact response (fuel removed) 
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(d) Calculated aircraft debris 
Figure 9-22. Calculated floor 96 contents, fuel, and aircraft debris distribution. 
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(c) Calculated impact response (fuel removed) 
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(d) Calculated aircraft debris 
Figure 9-23. Calculated floor 97 contents, fuel, and aircraft debris distribution. 
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(c) Calculated impact response (fuel removed) 



(d) Calculated aircraft debris 
Figure 9-24. Calculated floor 98 contents, fuel, and aircraft debris distribution. 
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The detailed predictions of the fuel dispersion and distribution using SPH in the global impact analyses 
had significant uncertainties in the absence of improved validation testing. However, some aspects of the 
distribution had a higher confidence. The floors confined the vertical motion of the fuel, and the floor-by- 
floor distribution of fuel was controlled more by the geometry of the tower and impact conditions. As a 
result, this distribution by floor has a higher level of confidence. Similarly, the interior contents and 
partition walls, and the damage to these structures, controlled the spread of fuel. 



Table 9-4. Fuel and aircraft debris distribution for the base case WTC 1 impact. 


Tower Location 


Aircraft Fuel 


Aircraft Debris 


Total Outside Tower 


6,700 lb 


17,400 lb 


WTC 1 Floor 92 


8101b 


260 1b 


WTC 1 Floor 93 


6,100 lb 


22,600 lb 


WTC 1 Floor 94 


16,100 lb 


96,000 lb 


WTC 1 Floor 95 


12,200 lb 


28,000 lb 


WTC 1 Floor 96 


11,700 lb 


19,400 lb 


WTC 1 Floor 97 


9,500 lb 


6,000 lb 


WTC 1 Floor 98 


2,200 lb 


6,000 lb 


WTC 1 Floor 99 


7701b 


901b 


Total Weight 


66,100 lb 


196,000 lb 



9.5 



WTC 2 BASE CASE IMPACT ANALYSIS 



This case is referred to as Case B for the remainder of the WTC Investigation reports. The WTC 2 base 
case impact analysis was performed for a 0.62 s duration following initial impact of the aircraft nose with 
the south exterior wall. A side view of the base case WTC 2 global impact model is shown in 
Figure 9-25. The tower interior contents were removed, and the tower structures were shown as 
transparent in the figure so that the impact response in the tower interior was visible. The floor numbers 
and core columns for the tower model are labeled in the figure. The corresponding side view is used to 
illustrate the base case global aircraft impact response in Figure 9-26. The response is shown at intervals 
of 0.1 s from impact through the initial 0.5 s of the response. Corresponding plan view images of the 
progression of aircraft impact response are shown in Figure 9-27. 
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Figure 9-25. WTC 2 base case global impact model. 
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(b) Time=0.10 s 
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(c) Time=0.20 s 
Figure 9-26. WTC 2 base case global impact analysis (side view). 
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(e) Time=0.40 s 
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(f) Time=0.50 s 
Figure 9-26. WTC 2 base case global impact analysis (side view) (continued). 
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(c) Time=0.20 s 
Figure 9-27. WTC 2 base case global impact analysis (plan view). 
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(d) Time=0.30 s 
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(e) Time=0.40 s 




(f) Time=0.50 s 
Figure 9-27. WTC 2 base case global impact analysis (plan view) (continued). 



NISTNCSTAR 1-2B, WTC Investigation 



233 



Chapter 9 



The initial 0.1 s of the base case global aircraft impact response, shown in Figure 9-26(b) and 
Figure 9-2 7(b), was dominated by the impact, penetration, and fragmentation of the forward fuselage 
structures. The engines and leading portions of the wings penetrated the exterior wall. The forward 
fuselage structures were severely damaged both from the penetration through the exterior columns and the 
penetration of the 81st floor slab that sliced the fuselage structures in half. The downward trajectory of 
the aircraft structures caused the airframe to collapse against the floor, and the subsequent debris motion 
was redirected inward along a more horizontal trajectory parallel to the floor. The higher impact speed 
and short truss floor span in this impact orientation had the forward fuselage structures well into the tower 
core by this time. 

At 0.2 s after impact, the full penetration of the aircraft into the tower was just completed, as shown in 
Figure 9-26(c) and Figure 9-27(c). The airframe was mostly broken up, but some large sections of the 
aft fuselage and tail were still intact, having penetrated through the opening in the tower wall produced by 
the forward fuselage structures. The aircraft fuel cloud was starting to spread out but was still relatively 
dense, and the leading edge of the fuel was approximately one-third through the tower core. By 0.2 s, the 
downward trajectory of the aircraft structures had transferred sufficient vertical load that the truss floor 
structures on the 80th and 81st floors were starting to collapse in the impact zone. 

At 0.3 s after impact, the aircraft fuel cloud had penetrated approximately two-thirds the distance through 
the core and was spreading out, as shown in Figure 9-26(d) and Figure 9-27(d). However, the 
subsequent motion of the aircraft fragments and fuel debris cloud began to be noticeably slowed beyond 
this time. The fuel and debris did continue to spread through the tower, but at a much slower rate, as seen 
in the remaining images in Figure 9-26 and Figure 9-27. The spread of the fuel and debris cloud was 
more rapid and extensive in the open truss floor regions than through the core as a result of the open 
volume above the workstations in the truss floor zone. 

The load transfer of the base case WTC 2 aircraft impact can be shown by the history of the resultant 
momentum as shown in Figure 9-28. The curve illustrates an initial rate of load transfer during the first 
0.1 s of impact as the forward fuselage penetrated the exterior wall and impacted the interior structures. 
Between 0.1 s and 0.2 s, a more rapid load transfer rate was observed, as the area of the impact was larger 
(extending outward in the wing impact regions), and a higher percentage of the aircraft mass was 
impacting the interior structures. At 0.2 s, as indicated in Figure 9-28, the aircraft completely penetrated 
the building and retained approximately 30 percent of its initial momentum. Beyond this time, the rate of 
load transfer was steadily decreasing with very little load transfer after approximately 0.4 s. The behavior 
was very similar to that of the base case WTC 1 impact, shown in Figure 9-6, but with a slightly 
compressed time scale resulting from the higher impact speed on WTC 2. 
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Figure 9-28. Normalized aircraft momentum for the WTC 2 base case impact. 

The aircraft was severely broken into debris as a result of the impact with WTC 2 as shown in 
Figure 9-29. At the end of the impact analysis, the aircraft was broken into thousands of debris fragments 
of various size and mass, as shown in Figure 9-29(b). A closer inspection of the debris field showed that 
larger fragments still occurred for specific components such as the engines and landing gear components. 
This behavior was very similar to the WTC 1 aircraft breakup shown in Figure 9-7. 



9.5.1 



WTC 2 Base Case Impact Analysis - Structural Damage 



The structural damage to the WTC 2 tower by the base case impact conditions is described in this section. 
The details of the base case WTC 2 aircraft impact conditions were summarized previously in Section 9.3 
(shown in Figure 9-2 and summarized in Table 9-2). The primary structural components of interest were 
the exterior wall, core columns and core framing components, and the floor structures and concrete floor 
slab. The response of the internal contents, walls, and partitions as well as fuel dispersion was also of 
interest as initial conditions for the fire analysis. 
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(a) Aircraft structure (time=0.00 s) 
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(b) Aircraft debris field (time=0.62 s) 
Figure 9-29. Aircraft breakup and debris in the base case WTC 2 impact. 
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Exterior Wall Damage 

A comparison of the south exterior wall observed and calculated damage from the base case WTC 2 
global impact analysis is shown in Figure 9-30. The calculated impact damage to the exterior wall is 
shown with color fringes representing plastic strain magnitude with undamaged sections in blue and 
strains at or above 5 percent shown in red. The schematic of observed damage was developed from 
inspections of the video and photographic data collected on the tower after impact. Both the observed and 
calculated damage regions shown in Figure 9-30 illustrate a region of the exterior wall from column 402 
to column 446, extending from floor 76 to floor 86 (spandrels at floors 77 through 86). 

The exterior wall completely failed in the regions of the fuselage, engine, and fuel-filled wing section 
impacts. Damage to the exterior wall extended to the wing tips, but the exterior columns were not 
completely failed in the outer wing and vertical stabilizer impact regions. Failure of the exterior columns 
occurred both at the bolted connections between column ends and at various locations in the column, 
depending on the local severity of the impact load and the proximity of the bolted connection to the 
impact. 

The initial observation from the comparison of the calculated and observed damage was that the geometry 
and location of the impact damage zone were in good agreement. This agreement in the position and 
shape of the impact damage served to validate the geometry of the aircraft model, including the aircraft 
orientation, trajectory, and flight distortions of the wings. The agreement of both the mode and 
magnitude of the impact damage served to partially validate the constitutive and damage modeling of the 
aircraft and exterior wall of the tower. 

Core Structural Damage 

The estimation of the damage to the core columns and core beams was important for the subsequent 
analysis of structural stability and collapse. The overall model for the WTC 2 core structure and 
calculated impact damage to the core is shown in Figure 9-31. The figure shows that the core had 
significant damage in the region close to the impact point, in particular the southeast corner of the core. 
The columns in line with the aircraft fuselage failed on the impact side, and several of the core beams 
were also severely damaged or failed in the impact zone. In some cases, failure of the column splices 
located on floors 77, 80, and 83 contributed significantly to the failure of the core columns. This was 
particularly true for the heavy column number 1001 at the corner of the core that failed at the three splice 
locations. 

The calculated damage to the core columns by row is shown in Figure 9-32. The columns are shown 
with color fringes representing plastic strain magnitude with undamaged sections in blue and strains at or 
above 5 percent shown in red. A summary of the column damage is listed in Table 9-5 and shown 
graphically in Figure 9-33. The qualitative classification of the column damage levels were provided 
previously in Figure 9-12. A summary of the floor-by-floor damage to the core framing is shown in 
Figure 9-34. 
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(b) Calculated damage (t=0.62 s) 
Figure 9-30. Base case impact damage to the WTC 2 exterior wall. 
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(a) Initial geometry 




(a) Calculated impact damage 
Figure 9-31. Base case impact damage to the WTC 2 core. 
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Figure 9-32. Base case impact damage to the WTC 2 core columns. 
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Table 9-5. Summary of core column damage for the base t 


:ase WTC 2 impact. 


Column 


Location 


Damage Level 


Lateral Deflection of 
Column Centerline (in.) 


Column 801 


Floor 79 


Heavy 


10 


Column 901 


Floors 79-82 


Severed 




Column 902 


Floor 79 


Heavy 


32 


Column 903 


Floors 77-83 


Severed 




Column 904 


Floor 79 


Moderate 




Column 905 


Floor 79 


Heavy 


18 


Column 1001 


Floors 77-83 


Severed 




Column 1002 


Floors 79-81 


Severed 




Column 1003 


Floor 80 


Severed 
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Floor 80 


Heavy 


18 
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Figure 9-33. Summary of core column damage for the base case WTC 2 impact. 
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(a) Floor 77 core framing damage 




(c) Floor 7S 
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(e) Floor 81 core framing damage 
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(b) Floor 78 core framing damage 




(d) Floor 80 core framing damage 
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(f) Floor 82 core framing damage 



Figure 9-34. Base case impact damage to the WTC 2 core floor framing (plan view). 
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Floor Truss and Slab Damage 

An overall view for the floor truss structure in the WTC 2 impact zone, along with the calculated base 
case impact damage to the trusses, is shown in Figure 9-35. The figure shows a front view of the trusses 
that were modeled in detail before and after impact by the aircraft. Note that the coarsely meshed trusses 
on these floors are not shown in the figure. The trusses experienced significant damage in the impact 
zone, with the largest damage on floor 81. A plan view of the calculated damage to the trusses on each 
floor is shown in Figure 9-36. The calculated impact response produced severe damage to the truss 
structures in the primary impact path of the fuselage. The truss structures were severely damaged from 
the exterior wall to the core. The truss floor system on floors 79 and 81 had sufficient damage from the 
impact that truss floor sections were sagging downward as a result of the impact. 




Floor SO 

Floor 79 

Floor 78 



(a) Initial detailed truss structures 



Floor 82 r^aSgrggSS 

_3 
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(b) Calculated damage 
Figure 9-35. Base case impact damage to the WTC 2 floor truss (front view). 

The calculated damage to the WTC 2 floor slab for floors 78 through 83 are shown in Figure 9-37. The 
fringes of damage were set such that the concrete failed in the regions colored red (2 percent plastic 
strain). In these regions, it is expected that the concrete had been severely damaged and potentially 
removed, exposing the supporting metal decking. The strength of the floor slab was severely reduced in 
the analysis beyond this strain to model the residual strength of the metal deck after the concrete failure, 
breakup, and removal. At a plastic strain of 30 percent, corresponding to failure levels for the metal 
decking material, the elements were eroded (seen as holes ruptured in the floor slab shown). The line of 
damage, visible around the edge of the core in Figure 9-37, was believed to be the result of numerical 
precision errors and interference in the modeling methodology at the connection of the floor slab rather 
than impact damage. 
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(a) Floor 79 truss damage 



(b) Floor 80 truss damage 
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(c) Floor 81 truss damage 




(d) Floor 82 truss damage 
Figure 9-36. Base case impact damage to the WTC 2 floor trusses (plan view). 
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(a) Floor 78 slab damage 
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(c) Floor 80 slab damage 
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(b) Floor 79 slab damage 



Impact 



^ ^ 



^p b r i u 



W 




mi 



rnrir 



__!_ J 



(d) Floor 81 slab damage 



(f) Floor 83 slab damage 



Figure 9-37. Base case impact damage to the WTC 2 floor slab (plan view). 
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9.5.2 WTC 2 Base Case Impact Analysis - Fuel and Debris Distributions 

The global impact results presented in this section include the distribution of the jet fuel and aircraft 
debris in the WTC 2 tower, and the damage to the building contents (partition walls and workstations). 
For the base case WTC 2 global impact analysis, the calculated distribution of the fuel in the tower and 
shape of the fuel cloud in a plan view and side view are shown in Figure 9-38 and Figure 9-39, 
respectively. At the end of the analysis, the residual momentum of the jet fuel in the impact direction was 
less than one percent of the initial momentum, indicating that the fuel cloud was nearly at rest at about 
0.62 s. 

To more clearly present the calculated response of the structures that influenced the fire propagation, the 
structural components were removed from the visualization, with the exception that the core columns 
were maintained in the visualizations for reference positions. A plan view of the response of the 
remaining building contents and aircraft debris are shown in a plan view in Figure 9-40. Similar plan 
views of floor-by-floor slices through the building contents and debris field are provided in Figure 9-41 
through Figure 9-46 for floors 78 through 83, respectively. The bulk of the aircraft debris and fuel was 
arrested prior to exiting the tower structures. However, a significant amount of aircraft debris was 
calculated to exit the north and east sides of the tower (Sides 300 and 200 of WTC 2). 

A quantitative characterization of the fuel and aircraft debris distribution was obtained by slicing the 
model at vertical floor locations and calculating the mass at each floor level. A summary of the floor-by- 
floor fuel and debris distributions is given in Table 9-6. The bulk of the fuel and aircraft debris was 
deposited in floors 78 through 80, with the greatest concentration of aircraft debris on floor 80, and the 
largest concentration of aircraft fuel on floors 79, 81, and 82. Approximately 14,000 lb, or 5 percent, of 
the total aircraft mass was eliminated from the debris cloud in the final state as a result of the erosion in 
the aircraft structures due to impact and breakup. This eroded mass was maintained in the calculation but 
eliminated from consideration in the contact algorithm. As a result, any residual momentum at the time 
of erosion could not be subsequently transferred to the tower. 

The calculated debris distribution included 55,800 lb of debris and 10,600 lb of aircraft fuel outside of the 
tower at the end of the impact analysis, either rebounding from the impact face or passing through the 
tower. These estimates of mass outside the tower were expected to be overestimated in the calculation 
since the exterior walls were not modeled with windows that could contain the fuel cloud and small debris 
inside the towers. In addition, the impact behavior of the aircraft fuel cloud did not include the ability to 
stick to, or wet, interior components. Rather, the aircraft fuel SPH particles tended to bounce off of 
internal structures (see Section 9.4.2). 
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(a) Calculated impact response (floor slab removed) 
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(b) SPH fuel distribution 
Figure 9-38. Calculated fuel distribution in the base case WTC 2 analysis (plan view). 
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(a) Calculated impact response 



Time = 0.621 9 
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(b) SPH fuel distribution 
Figure 9-39. Calculated fuel distribution in the base case WTC 2 analysis (side view). 
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(a) Impact configuration (plan view) 
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(b) Calculated impact response (t=0.622) 
Figure 9-40. Calculated WTC 2 building contents, fuel, and aircraft debris distribution. 
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(a) Initial impact configuration 
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(b) Calculated impact response 
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(c) Calculated impact response (fuel removed) 
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(d) Calculated aircraft debris 
Figure 9-41. Calculated floor 78 contents, fuel, and aircraft debris distribution. 
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(a) Initial impact configuration 
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(b) Calculated impact response 
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(c) Calculated impact response (fuel removed) 
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(d) Calculated aircraft debris 
Figure 9-42. Calculated floor 79 contents, fuel, and aircraft debris distribution. 
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(a) Initial impact configuration 



"!-.'..--..-lj- ?'■.." f 




(b) Calculated impact response 
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(c) Calculated impact response (fuel removed) 






(d) Calculated aircraft debris 
Figure 9-43. Calculated floor 80 contents, fuel, and aircraft debris distribution. 
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(a) Initial impact configuration 




(b) Calculated impact response 




(c) Calculated impact response (fuel removed) 
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(d) Calculated aircraft debris 
Figure 9-44. Calculated floor 81 contents, fuel, and aircraft debris distribution. 
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(b) Calculated impact response 
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(c) Calculated impact response (fuel removed) 



(d) Calculated aircraft debris 
Figure 9-45. Calculated floor 82 contents, fuel, and aircraft debris distribution. 
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(a) Initial impact configuration 
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(b) Calculated impact response 






(c) Calculated impact response (fuel removed) 



(d) Calculated aircraft debris 
Figure 9-46. Calculated floor 83 contents, fuel, and aircraft debris distribution. 
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Table 9-6. Fuel and aircraft debris distribution for the base case WTC 2 impact. 


Tower Location 


Aircraft Fuel 


Aircraft Debris 


Total Outside Tower 


10,600 lb 


55,800 lb 


WTC 2 Floor 77 


1,300 lb 


400 1b 


WTC 2 Floor 78 


6,200 lb 


4,800 lb 


WTC 2 Floor 79 


11,400 lb 


16,200 lb 


WTC 2 Floor 80 


6,000 lb 


83,800 lb 


WTC 2 Floor 81 


14,400 lb 


27,300 lb 


WTC 2 Floor 82 


10,600 lb 


3,600 lb 


WTC 2 Floor 83 


1,500 lb 


4,300 lb 


WTC 2 Floor 84 


2001b 


500 1b 


Total Weight 


62,000 lb 


197,600 lb 



9.6 



ADDITIONAL WTC 1 GLOBAL IMPACT ANALYSES 



In addition to the base case impact analyses, additional impact analyses were performed for each tower to 
provide a range of calculated impact-induced damage. The variations in impact analysis parameters were 
developed based on the result of the uncertainty analyses described in Chapter 8 and additional 
evaluations of the parameter uncertainties. Additional guidance on the combination of parameters was 
obtained by engineering analysis of expected responses. 

Two additional global impact analyses were performed for each WTC tower impact to evaluate a more 
severe and less severe impact scenario. The combinations of impact conditions for the more and less 
severe WTC 1 impact scenarios are compared to the corresponding parameters in the base case analysis in 
Table 9-7. The impact speeds were 414 mph and 472 mph in the less severe and more severe impact 
scenarios, respectively. These were the upper and lower bounds obtained from the analysis of aircraft 
impact conditions described in Chapter 7. The aircraft vertical trajectory angle was also varied from 
13.6 degrees to 7.6 degrees for the less severe and more severe impact scenarios, respectively. This 
resulted in more impact energy directed inward toward the core in the more severe impact scenario and 
less impact energy imparted to the core for the less severe impact. The lateral trajectory was not varied 
since the impact was close to centered on the tower and normal to the north face of WTC 1. A small 
variation in the lateral approach angle would have little effect on the energy of the aircraft debris entering 
the tower and core. 

The parameters varied for the aircraft model were the weight of the aircraft and the ductility of the aircraft 
materials. A 5 percent variation on the total aircraft weight was considered, with the additional weight 
applied for the more severe impact. The failure strain was varied from 75 to 125 percent of the baseline 
value for the less severe and more severe impact scenarios, respectively. This relatively large variation in 
aircraft material ductility was applied for multiple reasons. First, no material characterization testing of 
specimens cut from a Boeing 767 were performed as a part of this investigation. All of the material 
properties used for the aircraft were obtained from sources available in the open literature, as described in 
Section 2.5. Second, the variation in ductility was used as the single parameter in these additional 
analyses to evaluate the uncertainties in aircraft material energy absorption capacity. An increase in 
aircraft material strength would have a similar effect to an increase in material ductility for producing 
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increased impact damage to the towers. Finally, the material failure parameters were influenced by the 
resolution of the models in the impact analysis. This effect was described previously in Section 2.3.1, and 
the material models were developed for the resolution in the global impact analysis. However, the mesh 
refinement effects introduced an increased uncertainty on the failure strains in these analyses. 

Table 9-7. Input parameters for additional WTC 1 global impact analyses. 



Analysis Parameters 


Base Case 


More Severe 


Less Severe 


Flight 
Parameters 


Impact speed 


443 mph 


472 mph 


414 mph 


Trajectory - pitch 


10.6° 


7.6° 


13.6° 


Trajectory - yaw 


0.0° 


0.0° 


0.0° 


Orientation - pitch 


8.6° 


5.6° 


11.6° 


Orientation - yaw 


0.0° 


0.0° 


0.0° 


Aircraft 
Parameters 


Weight 


100 percent 


105 percent 


95 percent 


Failure Strain 


100 percent 


125 percent 


75 percent 


Tower 
Parameters 


Failure Strain 


100 percent 


80 percent 


120 percent 


Live Load Weight 3 


25 percent 


20 percent 


25 percent 



a. Live load weight expressed as a percentage of the design live load. 

Finally, the parameters varied for the tower model were the ductility of the steel used in the tower 
construction and the weight of the contents inside the tower. A variation of 20 percent was used to 
account for the uncertainty in failure strain for the tower materials. However, for the tower the material, 
failure strain was increased for the less severe impact analysis. The combination of reducing all of the 
aircraft material ductilities by 25 percent and increasing the tower material ductilites by 20 percent for the 
less severe impact scenario covered a wide range in relative aircraft and tower strength assumptions. 

The variations in internal tower contents (live load weight in Table 9-7) are specified as a percentage of 
the design live load as described in Section 3.2.4. 



9.6.1 



More Severe WTC 1 Global Impact Analysis 



This case is referred to as Case C for the remainder of the WTC investigation reports. A side view of the 
more severe WTC 1 global impact response is shown in Figure 9-47. The tower interior contents were 
removed and the tower structures were shown as transparent so that the impact response in the tower 
interior can be seen. Similar plan view images of the progression of aircraft impact response are shown in 
Figure 9—48. The corresponding views were used previously to illustrate the base case global aircraft 
impact response in Figure 9-4 and Figure 9-5, respectively. 

The more severe impact response is qualitatively very similar to that of the base case impact analysis. 
The initial 0.1 s of the base case global aircraft impact response, shown in Figure 9-47(b) and 
Figure 9^18(b), was dominated by the impact, penetration, and fragmentation of the forward fuselage 
structures. The engines and forward wing sections have impacted the exterior wall. At 0.2 s after impact, 
the wings had completely penetrated the exterior wall and only the tail structures were still outside of the 
tower, as shown in Figure 9-47(c) and Figure 9-48(c). The wing structures were completely fragmented 
by the penetration through the exterior wall. The aircraft fuel cloud was starting to spread out, but was 
still relatively dense and the leading edge of the fuel was entering the tower core. At 0.3 s after impact, 
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the aircraft was completely inside of the tower (full penetration completed at approximately 0.24 s), as 
shown in Figure 9-47(d) and Figure 9-48(d). The aircraft fuel cloud penetrated approximately three- 
quarters of the distance through the core and was spreading out. However, the subsequent motion of the 
aircraft fragments and fuel debris cloud began to be noticeably slowed beyond this time. The fuel and 
debris did continue to spread through the tower but at a much slower rate as seen in the remaining images 
in Figure 9-47 and Figure 9-48. 

Exterior Wall Damage for the More Severe WTC 1 Impact 

A comparison of the north exterior wall observed and calculated damage from the more severe WTC 1 
global impact analysis is shown in Figure 9-49. The calculated impact damage to the exterior wall is 
shown with color fringes representing plastic strain magnitude, with undamaged sections in blue and 
strains at or above 5 percent shown in red. The wall regions shown in Figure 9^19 illustrate a region of 
the exterior wall from column 108 to column 152, extending from floor 91 to floor 100 (spandrels at 
floors 92 through 100). The calculated and observed magnitude and mode of impact damage on the 
exterior wall was still in good agreement for this more severe impact analysis. 

Core Structural Damage for the More Severe WTC 1 Impact 

The overall model for the WTC 1 core structure and calculated response for the more severe global 
impact analysis is shown in Figure 9-50. The figure shows that the core had extensive damage in the 
region close to the impact point. The columns in line with the aircraft fuselage failed on the impact side, 
and several of the core beams were also severely damaged or failed in the impact zone. In some cases, 
failure of the column splices located on floors 92, 95, and 98 contributed significantly to the failure of the 
core columns. 

The calculated damage to the core columns by row is shown in Figure 9-51, and the floor-by-floor 
damage to the core framing is shown in Figure 9-52. The core structures are shown with color fringes 
representing plastic strain magnitude, with undamaged sections in blue and strains at or above 5 percent 
shown in red. A summary of the core column damage is provided in Table 9-8 and shown graphically in 
Figure 9-53, with the qualitative classification of the column damage levels provided previously in 
Figure 9-12. A total of six columns were severed and three columns were heavily damaged, compared to 
three columns severed and four columns heavily damaged in the base case WTC 1 impact analysis. 
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Figure 9-47. More severe WTC 1 global impact analysis (side view). 
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(f) Time=0.50 s 
Figure 9-47. More severe WTC 1 global impact analysis (side view) (continued). 
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(c) Time=0.20 s 
Figure 9-48. More severe WTC 1 global impact analysis (plan view). 
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(f) Time=0.50 s 
Figure 9-48. More severe WTC 1 global impact analysis (plan view) (continued). 
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(b) Calculated damage (t=0.685 s) 
Figure 9-49. More severe impact damage to the WTC 1 exterior wall. 
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(a) Initial geometry 




(a) Calculated impact damage 
Figure 9-50. More severe impact response of the WTC 1 core. 
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(a) Columns 501-1001 




(d) Columns 504-1004 
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(e) Columns 505-1005 
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Figure 9-51. More severe impact response of the WTC 1 core columns. 
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(a) Floor 93 core framing damage 




(c) Floor 95 core framing damage 
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(b) Floor 94 core framing damage 




(d) Floor 96 core framing damage 
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(e) Floor 97 core framing damage (f) Floor 98 core framing damage 

Figure 9-52. More severe impact damage to the WTC 1 core floor framing (plan view). 
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Table 9-8. Core column response for the more severe WTC 1 impact. 


Column 


Location 


Damage Level 


Lateral Deflection of 
Column Centerline (in.) 


Column 503 


Floor 95-96 


Severed 




Column 504 


Floors 92-96 


Severed 




Column 505 


Floors 93-96 


Severed 




Column 506 


Floors 93-95 


Heavy 


24 


Column 603 


Floors 96-97 


Moderate 




Column 604 


Floors 92-96 


Severed 




Column 605 


Floors 94-95 


Moderate 




Column 606 


Floors 94 


Light 




Column 702 


Floor 97 


Light 




Column 703 


Floor 96 


Moderate 




Column 704 


Floors 92-96 


Severed 




Column 705 


Floor 95 


Moderate 




Column 706 


Floors 93-95 


Severed 




Column 802 


Floor 96 


Light 




Column 803 


Floors 96-97 


Moderate 




Column 804 


Floor 94-96 


Moderate 




Column 805 


Floors 93-95 


Heavy 


20 


Column 903 


Floor 96 


Light 




Column 904 


Floors 95-96 


Heavy 


19 


Column 905 


Floor 95 


Light 
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Figure 9-53. Summary of core column damage for the more severe WTC 1 impact. 

Floor Truss and Slab Damage for the More Severe WTC 1 Impact 

An overall view of the floor truss structure in the impact zone, along with the calculated more severe 
impact damage to the floor trusses, is shown in Figure 9-54. The figure shows that the trusses 
experienced significant damage in the impact zone. A plan view of the calculated damage to the truss on 
each floor is shown in Figure 9-55. The calculated impact response produced severe damage to the truss 
structures in the primary impact path of the fuselage. The truss structures were severely damaged from 
the exterior wall to the core. The truss floor system on floors 94 through 96 were damaged and sagging 
downward as a result of the impact loading. 

The magnitude of truss floor damage was very similar for the base case and more severe global impact 
analyses. When the floor-by-floor damage was compared for the base case and more severe impact 
analyses, Figure 9-15 and Figure 9-55 respectively, the damage appears to be slightly less for the more 
severe impact analysis. The parameters used in the more severe global impact analysis would primarily 
contribute to an increased damage magnitude for the tower structures. However, the downward impact 
trajectory angle was reduced from the 10.6 degree angle in the base case analysis to a 7.6 degree angle in 
the more severe impact analysis. This would have the effect of directing more of the impact energy 
inward toward the tower core but reducing the normal downward force on the floor structures in the 
impact zone. As a result, the combined effects of the analysis parameter variations produced slightly less 
damage to the truss structure in the more severe impact analysis scenario. 
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(b) Calculated damage (t=0.685 s) 
Figure 9-54. More severe impact damage to the WTC 1 floor truss (front view). 
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(c) Floor 96 truss damage (d) Floor 97 truss damage 

Figure 9-55. More severe impact damage to the WTC 1 floor truss (plan view). 
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The calculated more severe impact damage to the WTC 1 floor slab for floors 94 through 97 are shown in 
Figure 9-56. The fringes of damage were set such that the concrete is failed in the regions colored red 
(2 percent plastic strain was used corresponding to the zero strength strain limit for the concrete in 
unconfined compression). At these strain levels, the concrete was severely damaged and probably 
removed, exposing the supporting metal decking. Beyond 2 percent plastic strain, the strength of the 
floor slab was severely reduced in the analyses to model the residual strength of the metal deck after the 
concrete failure, breakup, and removal. At a plastic strain of 30 percent, corresponding to failure levels 
for the metal decking material, the elements were eroded (seen as holes ruptured in the floor slab shown). 
The line of damage, visible around the edge of the core, was believed to be the result of numerical 
precision errors and interference in the modeling methodology at the connection of the floor slab rather 
than impact damage. 
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Figure 9-56. More severe impact damage to the WTC 1 floor slabs (plan view). 
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The magnitude of floor slab damage was, in general, very similar for the base case and more severe global 
impact analyses. When the floor-by-floor damage was compared for the base case and more severe 
impact analyses, Figure 9-16 and Figure 9-56 respectively, the damage appears to be slightly less for the 
more severe impact analysis. Similar to the truss damage, the reduced damage in the floor slab was 
believed to be the result of the reduction in the downward impact trajectory angle from 10.6 to 
7.6 degrees in the more severe impact analysis, reducing the normal downward force on the floor 
structures. 

Damage to WTC 1 Contents for the More Severe Impact 

The calculated damage to the WTC 1 contents for the more severe impact is shown in plan views in 
Figure 9-57 through Figure 9-61 for floors 94 through 98, respectively. A comparison to the calculated 
damage for the base case WTC 1 impact analysis, shown previously in Figure 9-20 through Figure 9-24, 
indicated that the content damage zone is very similar in width but extended further south through the 
tower in the more severe impact. The more severe impact produced significantly greater content damage 
on the far side of the core and extended more fully through the tower. 
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(b) Calculated response of contents 
Figure 9-57. Calculated more severe WTC 1 impact response of floor 94 contents. 
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(b) Calculated response of contents 
Figure 9-58. Calculated more severe WTC 1 impact response of floor 95 contents. 
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(b) Calculated response of contents 
Figure 9-59. Calculated more severe WTC 1 impact response of floor 96 contents. 
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(b) Calculated response of contents 
Figure 9-60. Calculated more severe WTC 1 impact response of floor 97 contents. 
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(b) Calculated response of contents 
Figure 9-61. Calculated more severe WTC 1 impact response of floor 98 contents. 

9.6.2 Less Severe WTC 1 Global Impact Analysis 

Exterior Wall Damage for the Less Severe WTC 1 Impact 

A comparison of the north exterior wall observed and calculated damage from the less severe WTC 1 
global impact analysis is shown in Figure 9-62. The calculated impact damage to the exterior wall is 
shown with color fringes representing plastic strain magnitude, with undamaged sections in blue and 
strains at or above 5 percent shown in red. The wall regions shown in Figure 9-62 are from column 108 
to column 152, extending from floor 91 to floor 100 (spandrels at floors 92 through 100). The magnitude 
and mode of impact damage on the exterior wall were still in good agreement with the observed damage 
for this less severe impact scenario. 
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(b) Calculated damage (t=0.735 s) 
Figure 9-62. Less severe impact damage to the WTC 1 exterior wall. 
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Core Structural Damage for the Less Severe WTC 1 Impact 

The overall model for the WTC 1 core structure and calculated response for the less severe global impact 
analysis is shown in Figure 9-63. The figure shows that the core had a limited damage confined to the 
region nearest to the impact point. The calculated damage to the core columns by row is shown in 
Figure 9-64, and the floor-by-floor damage to the core framing is shown in Figure 9-65. The core 
structures are shown with color fringes representing plastic strain magnitude, with undamaged sections in 
blue and strains at or above 5 percent shown in red. A summary of the column damage is provided in 
Table 9-9 and shown graphically in Figure 9-66, with the qualitative classification of the column damage 
levels provided previously in Figure 9-12. Only one column was severed, and two columns were heavily 
damaged, compared to three severed columns and four heavily damaged columns in the base case WTC 1 
impact analysis. The failure of the column splices located on floors 92 and 95 contributed to the failure of 
the core column. 



Table 9-9. Core column response 


for the less severe Wl 


rc 1 impact. 


Column 


Location 


Damage Level 


Lateral Deflection of 
Column Centerline (in.) 


Column 504 


Floors 92-96 


Severed 




Column 505 


Floors 93-96 


Heavy 


14 


Column 506 


Floors 93-94 


Moderate 




Column 604 


Floor 96 


Light 




Column 606 


Floors 94 


Light 




Column 702 


Floor 96 


Moderate 




Column 703 


Floor 96 


Light 




Column 704 


Floor 96 


Heavy 


12 


Column 705 


Floor 96 


Light 




Column 802 


Floor 96 


Moderate 




Column 805 


Floor 96 


Light 
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(a) Initial geometry 




(a) Calculated impact damage 
Figure 9-63. Less severe impact response of the WTC 1 core. 
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Figure 9-64. Less severe impact response of the WTC 1 core columns. 
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(e) Floor 97 core framing damage (f) Floor 98 core framing damage 

Figure 9-65. Less severe impact damage to the WTC 1 core floor framing (plan view). 



NISTNCSTAR 1-2B, WTC Investigation 



281 



Chapter 9 




4 



M 



508 



1008 



© 



0© 



501 



1001 



Severed •Heavy OModerate 



Figure 9-66. Summary of core column damage for the less severe WTC 1 impact. 

Floor Truss and Slab Damage for the Less Severe WTC 1 Impact 

An overall view of the floor truss structure in the impact zone, along with the calculated less severe 
impact damage to the floor trusses, is shown in Figure 9-67. The figure shows that the trusses 
experienced significant damage in the impact zone. A plan view of the calculated damage to the truss on 
each floor is shown in Figure 9-68. The calculated impact response produced severe damage to the truss 
structures in the primary impact path of the fuselage. The truss structures were severely damaged from 
the exterior wall to the core. The truss floor system on floors 94 through 96 were damaged and sagged 
downward as a result of the impact loading. 
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(b) Calculated damage (t=0.735 s) 
Figure 9-67. Less severe impact damage to the WTC 1 floor truss (front view). 

The magnitude of truss floor damage was very similar for the base case and more severe global impact 
analyses. When the floor-by-floor damage was compared for the base case and less severe impact 
analyses, Figure 9-15 and Figure 9-68, respectively, the damage was slightly increased for the less severe 
impact analysis. The parameters used in the less severe global impact analysis would primarily contribute 
to a reduced damage magnitude for the tower structures. However, the downward impact trajectory angle 
was increased from the 10.6 degree angle in the base case analysis to a 13.6 degree angle in the more 
severe impact analysis. This would have the effect of directing more of the impact energy downward, 
increasing the normal force on the floor structures in the impact zone. As a result, the combined effects 
of the analysis parameter variations produced a small increase in the damage to the truss structure in the 
less severe impact analysis scenario. 
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(c) Floor 96 truss damage (d) Floor 97 truss damage 

Figure 9-68. Less severe impact damage to the WTC 1 floor truss (plan view). 

The calculated less severe impact damage to the WTC 1 floor slab for floors 94 through 97 are shown in 
Figure 9-69. The magnitude of floor slab damage was very similar for the base case and less severe 
global impact analyses. When the floor-by-floor damage was compared for the base case and less severe 
impact analyses, Figure 9-16 and Figure 9-69, respectively, the damage appeared to be slightly increased 
in the truss floor impact zone for the less severe impact. Similar to the truss damage, the increased 
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damage in the floor slab was believed to be the result of the increase in the downward impact trajectory 
angle from 10.6 to 13.6 degrees in the more severe impact analysis, increasing the normal downward 
force on the floor structures. In contrast to the truss floor slab, the core floor slab damage was reduced in 
the less severe impact analysis. This was a result of less impact energy imparted to the core frame. In the 
base case analysis, the core frame was more severely damaged, resulting in coupled damage to the core 
slab. 
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(c) Floor 96 slab damage (d) Floor 97 slab damage 

Figure 9-69. Less severe impact damage to the WTC 1 floor slab (plan view). 

Damage to WTC 1 Contents for the Less Severe Impact 

The calculated damage to the WTC 1 contents for the less severe impact is shown in plan views in 
Figure 9-70 through Figure 9-74 for floors 94 through 98, respectively. A comparison to the calculated 
damage for the base case WTC 1 impact analysis, shown previously in Figure 9-20 through Figure 9-24, 
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indicated that the building contents damage zone was very similar in width but did not extend as far 
through the tower in the less severe impact. The less severe impact produced little content damage on the 
far side of the core and did not extend fully through the tower. Little or no debris penetration of the south 
wall of the tower was expected for the less severe impact condition. 
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(b) Calculated response of contents 
Figure 9-70. Calculated less severe WTC 1 impact response of floor 94 contents. 
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(a) Initial impact configuration 
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(b) Calculated response of contents 
Figure 9-71. Calculated less severe WTC 1 impact response of floor 95 contents. 
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(a) Initial impact configuration 




(b) Calculated response of contents 
Figure 9-72. Calculated less severe WTC 1 impact response of floor 96 contents. 



288 



NISTNCSTAR 1-2B, WTC Investigation 



Global Impact Analyses 




i m n m n 

JUJUUJU 

imnmn 



E 3 



in 



ru 



oh 


L 






f\ 


1 . 




L. 




X- J 

,. 1 


V 




, 


J , 




, . 



nmnmn 

UIUIJILTU 




nmnmn 



[Fufjij UTUTJIU^ 



(a) Initial impact configuration 
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(b) Calculated response of contents 
Figure 9-73. Calculated less severe WTC 1 impact response of floor 97 contents. 
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(b) Calculated response of contents 
Figure 9-74. Calculated less severe WTC 1 impact response of floor 98 contents. 

ADDITIONAL WTC 2 GLOBAL IMPACT ANALYSES 



9.7 



Two additional global impact analyses were performed for the WTC 2 impact to evaluate a more severe 
and less severe impact scenario. The combinations of impact conditions for the more and less severe 
impact scenarios are compared to the corresponding parameters in the base case analysis in Table 9-10. 
The impact speeds were 521 mph and 570 mph in the less severe and more severe impact scenarios, 
respectively. These were the upper and lower bounds obtained from the analysis of aircraft impact 
conditions described in Chapter 7. The aircraft vertical trajectory angle was also varied from 8 to 
5 degrees for the less severe and more severe impact scenarios, respectively. This resulted in more 
impact energy directed inward toward the core in the more severe impact scenario and less impact energy 
in the core for the less severe impact. The lateral trajectory was not varied in these analyses so that the 
starboard engine trajectory was aligned with exiting the northeast corner of the tower, as was observed 
from photographic evidence (see Section 9.11). 
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Table 9-10. Input parameters for additional WTC 2 global impact analyses. 



Analysis Parameters 


Base Case 


More Severe 


Less Severe 


Flight 
Parameters 


Impact Velocity 


546 mph 


570 mph 


521 mph 


Trajectory - pitch 


6.0° 


5.0° 


8.0° 


Trajectory - yaw 


13.0° 


13.0° 


13.0° 


Orientation - pitch 


5.0° 


4.0° 


7.0° 


Orientation - yaw 


10.0° 


10.0° 


10.0° 


Aircraft 
Parameters 


Weight 


100 percent 


105 percent 


95 percent 


Failure Strain 


100 percent 


115 percent 


75 percent 


Tower 
Parameters 


Contents Strength 


100 percent 


80 percent 


100 percent 


Failure Strain 


100 percent 


90 percent 


120 percent 


Live Load Weight" 


25 percent 


20 percent 


25 percent 



a. Live load weight expressed as a percentage of the design live load. 

The parameters varied for the aircraft model were the weight of the aircraft and the ductility of the aircraft 
materials. A 5 percent variation of the total aircraft weight was considered, with the additional weight 
applied for the more severe impact. The failure strain was varied from 75 to 115 percent of the baseline 
value for the less severe and more severe impact scenarios, respectively. This variation in aircraft 
material ductility was used as the single parameter in these additional analyses to evaluate the 
uncertainties in aircraft material energy absorption capacity. 

The parameters varied for the tower model were the ductility of the steel used in the tower construction 
and the weight and strength of the contents inside the tower. The less severe impact analysis used a 
20 percent increase in failure strain for the tower materials and the baseline values for the contents weight 
and strength. The more severe impact analysis used a 10 percent reduction in tower material ductility 
with a corresponding 20 percent reduction in strength of the internal contents and a reduction of the 
weight of the content from 25 percent to 20 percent of the design live load capacity for the truss and core 
floor areas. The design live load is described in Section 3.2.4. 

For WTC 2, the variations in the parameters from the base case were similar to those for WTC 1 (see 
Table 9-7), with two exceptions. The first exception was introduction of the strength of the building 
contents as a parameter. There was less information available about the layout of building contents in the 
WTC 2 impact zone and therefore a larger uncertainty associated with the contents was assumed. Recall 
that the workstation layout from WTC 1 was used for WTC 2. Thus, in the more severe case, the contents 
strength was reduced to 80 percent of the baseline value. 

The second exception was the failure strains for the aircraft and tower materials. For the more severe 
WTC 1 analysis, 125 percent and 80 percent of the baseline values were used for the aircraft and tower 
failure strains, respectively. For the more severe WTC 2 analysis, 115 percent and 90 percent of the 
baseline values were used. The more severe WTC 2 analysis was the final global impact analysis 
performed. Based on the previous analyses, the variation in damage levels indicated that the WTC 2 
more severe impact analysis would produce impact damage state that was not viable (e.g., the amount of 
debris exiting the north wall). To ensure that a viable damage state was obtained, the aircraft and tower 
materials were adjusted to the values presented in Table 9-10. 
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9.7.1 More Severe WTC 2 Global Impact Analyses 

This case is referred to as Case D for the remainder of the WTC investigation reports. A side view of the 
more severe WTC 2 global impact model is shown in Figure 9-75. The tower interior contents were 
removed and the tower structures were shown as transparent in the figure so that the impact response in 
the tower interior was visible. Similar plan view images of the progression of aircraft impact response are 
shown in Figure 9-76. The response is shown at intervals of 0.1 s from impact through the initial 0.5 s of 
the response. Corresponding images of the base case WTC 2 impact response were shown previously in 
Figure 9-26 and Figure 9-27, respectively. 

The initial 0.1 s of the more severe global aircraft impact response, shown in Figure 9-75(b) and 
Figure 9-76(b), was dominated by the impact, penetration, and fragmentation of the forward fuselage 
structures. The engines and forward wing structures also penetrated the exterior wall. At 0.2 s after 
impact, the full penetration of the aircraft into the tower was completed, as shown in Figure 9-75(c) and 
Figure 9-76(c). The airframe was mostly broken up, but some large sections of the aft fuselage and tail 
were still intact, having penetrated through the opening in the tower wall produced by the forward 
fuselage structures. The aircraft fuel cloud was starting to spread out, but was still relatively dense and 
the leading edge of the fuel was approximately one-half through the tower core. At 0.3 s after impact, the 
aircraft fuel cloud had penetrated the distance through the core and was spreading out, as shown in 
Figure 9-75(d) and Figure 9-76(d). The subsequent motion of the aircraft fragments and fuel debris 
cloud began to be noticeably slowed beyond this time. The slower spread of fuel and debris through the 
tower is seen in the remaining images in Figure 9-75 and Figure 9-76. The spread of the fuel and debris 
cloud was more rapid and extensive in the open truss floor regions than through the core as a result of the 
open volume above the workstations in the truss floor zone. 

Exterior Wall Damage for the More Severe WTC 2 Impact 

A comparison of the south exterior wall observed and calculated damage from the more severe WTC 2 
global impact analysis is shown in Figure 9-77. The calculated impact damage to the exterior wall is 
shown with color fringes representing plastic strain magnitude, with undamaged sections in blue and 
strains at or above 5 percent shown in red. The wall region shown in Figure 9-77 includes column 402 to 
column 446, extending from floor 76 to floor 86 (spandrels at floors 77 through 86). The mode and 
magnitude of the calculated and observed impact damage on the exterior wall were still in good 
agreement in this more severe impact analysis. 
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(a) Time=0.00 s 




(b) Time=0.10 s 




(c) Time=0.20 s 
Figure 9-75. More severe WTC 2 global impact analysis (side view). 
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(d) Time=0.30 s 
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(e) Time=0.40 s 
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(f) Time=0.50 s 
Figure 9-75. More severe WTC 2 global impact analysis (side view) (continued). 
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(b) Time=0.10 s 
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(c) Time=0.20 s 
Figure 9-76. More severe WTC 2 global impact analysis (plan view). 
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(f) Time=0.50 s 
Figure 9-76. More severe WTC 2 global impact analysis (plan view) (continued). 
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(a) Schematic of observed damage 
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Figure 9-77. More severe impact damage to the WTC 2 exterior wall. 

Core Structural Damage for the More Severe WTC 2 Impact 

The overall model for the WTC 2 core structure and calculated response for the more severe damage 
analysis is shown in Figure 9-78. The figure shows that the core had extensive damage in the region 
close to the impact point. The columns in line with the aircraft fuselage failed on the impact side, and 
several of the core beams were also severely damaged or failed in the impact zone. In some cases, failure 
of the column splices located on floors 77, 80, and 83 contributed significantly to the failure of the core 
columns. 
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The calculated damage to the core columns by row is shown in Figure 9-79, and the floor-by-floor 
damage to the core framing is shown in Figure 9-80. The core structures are shown with color fringes 
representing plastic strain magnitude, with undamaged sections in blue and strains at or above 5 percent 
shown in red. A summary of the column damage is provided in Table 9-11 and shown graphically in 
Figure 9-81, with the qualitative classification of the column damage levels provided previously in 
Figure 9-12. A total of 10 columns were severed, and one column was heavily damaged, compared to 
five columns severed and four columns heavily damaged in the base case WTC 2 impact analysis. 




(a) Initial geometry 




(a) Calculated impact damage 
Figure 9-78. More severe impact response of the WTC 2 core. 
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(a) Columns 1001-1008 





(b) Columns 901-908 




(c) Columns 801-807 (d) Columns 701-708 

Figure 9-79. More severe impact response of the WTC 2 core columns. 
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(a) Floor 78 core framing damage 




(c) Floor 80 core framing damage 
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(e) Floor 82 core framing damage 
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(b) Floor 79 core framing damage 






(d) Floor 81 core framing damage 
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(f) Floor 83 core framing damage 



Figure 9-80. More severe impact damage to the WTC 2 core floor framing (plan view). 
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Table 9-11. Core column response 


for the more severe WTC 2 impact. 


Column 


Location 


Damage Level 


Lateral Deflection of 
Column Centerline (in.) 


Column 602 


Floor 79 


Moderate 




Column 605 


Floor 79 


Moderate 




Column 701 


Floors 79-80 


Severed 




Column 702 


Floor 79 


Heavy 


16 


Column 703 


Floor 79 


Moderate 




Column 704 


Floor 79 


Light 




Column 705 


Floors 78-79 


Light 




Column 705 


Floor 78 


Light 




Column 801 


Floors 79-80 


Severed 




Column 802 


Floors 77-80 


Severed 




Column 803 


Floors 77-80 


Severed 




Column 804 


Floor 79 


Light 




Column 901 


Floors 80-81 


Severed 




Column 902 


Floor 79 


Moderate 




Column 903 


Floors 77-83 


Severed 




Column 904 


Floors 79-81 


Moderate 




Column 905 


Floors 79 & 81 


Light 




Column 907 


Floor 81 


Light 




Column 1001 


Floors 77-83 


Severed 




Column 1002 


Floors 79-83 


Severed 




Column 1003 


Floors 79-83 


Severed 




Column 1004 


Floors 79-83 


Severed 




Column 1005 


Floors 79-81 


Moderate 
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Figure 9-81. Summary of core column damage for the more severe WTC 2 impact. 

Strain Rate Observations for the More Severe WTC 2 Impact 

The more severe WTC 2 impact analysis had the highest impact speed of any of the global impact 
analyses. As a result, the more severe WTC 2 impact is an appropriate case for bounding the expected 
strain rates for the deformation of the WTC tower structural materials. Fringes of Von Mises effective 
strain rate for the south exterior wall and core framing are shown in Figure 9-82 and Figure 9-83, 
respectively. The corresponding times shown are at 90 ms and 200 ms for the exterior wall and core 
framing, respectively. These times were selected as representative of instants when rapid impact 
deformations over a significant portion of the structure were observed in the calculation. 

The strain rates in the exterior wall were significantly higher than those in the core. Peak strain rates were 
calculated at specific locations for short durations. These strain rates were on the order of a few thousand 
per second. However, these high strain rates were localized to small regions of the structure and for very 
short durations. Figure 9-82 shows that the majority of the deformations occurred with strain rates at or 
below approximately 100 per second. Similarly, the localized peak strain rates in the core were a few 
hundred per second, and Figure 9-83 shows that the majority of the deformation in the core occurred with 
strain rates at or below approximately 20 per second. 

The expected influence of these elevated strain rates on strength is small. The dominant strain rates seen 
in the figures are between 1 and 100 per second. The high rate material testing, described previously in 
Section 2.2.4, shows that the yield stress for the WTC tower materials typically increased by less than 
10 percent at these calculated rates. 
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Figure 9-82. Representative exterior wall strain rates for the more severe WTC 2 impact 

(view from tower exterior). 
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Figure 9-83. Representative core frame strain rates for the more severe WTC 2 impact 

(impact from the left of the figure). 
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Floor Truss and Slab Damage for the More Severe WTC 2 Impact 

An overall view of the floor truss structure in the WTC 2 impact zone, along with the calculated more 
severe impact damage to the trusses, is shown in Figure 9-84. The figure shows that the trusses 
experienced significant damage in the impact zone, with the heaviest damage on floor 81. A plan view of 
the calculated damage to the trusses on each floor is shown in Figure 9-85. The calculated impact 
response produced severe damage to the truss structures in the primary impact path of the fuselage. The 
truss structures were severely damaged from the exterior wall to the core. The truss floor system on 
floors 79 through 82 had sufficient damage from the impact that portions of the truss floor sections 
sagged downward as a result of the impact. 



Floor S2 



Floor SO 
Floor 79 

Floor 78 



(a) Initial detailed truss structures 




(b) Calculated damage (t=0.58 s) 
Figure 9-84. More severe impact damage to the WTC 2 floor truss (front view). 

The magnitude of truss floor damage was very similar for the base case and more severe global impact 
analyses. The parameters used in the more severe global impact analysis would primarily contribute to 
increased damage for the tower structures. However, the downward impact trajectory angle was reduced 
from the 6 degree angle in the base case analysis to a 5 degree angle in the more severe impact analysis. 
This would have the effect of directing more of the impact energy inward toward the tower core, but 
reducing the normal downward force on the floor structures in the impact zone. As a result, the combined 
effects of the analysis parameter variations produced very similar damage to the truss structure. 

The calculated damage to the WTC 2 floor slab for floors 78 through 83 for the more severe impact is 
shown in Figure 9-86. The fringes of damage were set such that the concrete failed in the regions colored 
red (2 percent plastic strain). In these regions, it is expected that the concrete slab was severely damaged 
and potentially removed, exposing the supporting metal decking. At a plastic strain of 30 percent, 
corresponding to failure levels for the metal decking material, the elements were eroded (seen as holes 
ruptured in the floor slab shown). The line of damage, visible around the edge of the core, was believed 
to be the result of numerical precision errors and interference in the modeling methodology at the 
connection of the floor slab rather than impact damage. 
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(b) Floor 80 truss damage 
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(d) Floor 82 truss damage 
Figure 9-85. More severe impact damage to the WTC 2 floor truss (plan view). 
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(a) Floor 78 slab damage 
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(c) Floor 80 slab damage 
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(d) Floor 81 slab damage 
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(e) Floor 82 slab damage (f) Floor 83 slab damage 

Figure 9-86. More severe impact damage to the WTC 2 floor slab (plan view). 
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Similar to the damage levels in the floor truss structures, the magnitude of floor slab damage was very 
similar for the base case and more severe global impact analyses. The parameters used in the more severe 
global impact analysis would primarily contribute to an increased damage magnitude for the tower 
structures. However, the downward impact trajectory angle was reduced from the 6 degree angle in the 
base case analysis to a 5 degree angle in the more severe impact analysis. This would have the effect of 
directing more impact energy inward toward the tower core, but reducing the normal downward force on 
the floor structures in the impact zone. As a result, the combined effects of the analysis parameter 
variations produced very similar damage magnitudes to the floor slab. 

Damage to WTC 2 Contents for the More Severe Impact 

The calculated damage to the WTC 2 contents for the more severe impact is shown in plan views in 
Figure 9-87 through Figure 9-92 for floors 78 through 83, respectively. A comparison to the calculated 
damage for the base case WTC 2 impact analysis, shown previously in Figure 9-41 through Figure 9-46, 
indicated that the tower contents damage zone was similar, with a slight increase in damage for the more 
severe impact. 
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(a) Initial impact configuration 



fhLX^> tJ i , nm; 

HLUUt> - ** ■ '^5 l^-K — = Ll LU LJ LU L 

rumnmn, 

LIILTlTLULI 

e^ run n m n nin n mnm q 

j U LU LJ UJ LJ LUTI ULULJIUU 

(b) Calculated response of contents 
Figure 9-87. Calculated more severe WTC 2 impact response of floor 78 contents. 
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(b) Calculated response of contents 
Figure 9-88. Calculated more severe WTC 2 impact response of floor 79 contents. 



308 



NISTNCSTAR 1-2B, WTC Investigation 



Global Impact Analyses 






* 




E 




p: : 




x 


i m n m r 

J 111 LJ JJJ L 

imrimr 






1 










[ 


3 






i 




i 


_£' I h 


h 








runn m n m n 

LRUUUTUUJU 


!H 


rLTTLTirrL 

LrUTLTLLT 


n^\l^ 


\\MJ^ 



lM3 IfE 



p El' 



(a) Initial impact configuration 
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(b) Calculated response of contents 
Figure 9-89. Calculated more severe WTC 2 impact response of floor 80 contents. 
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(a) Initial impact configuration 
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(b) Calculated response of contents 
Figure 9-90. Calculated more severe WTC 2 impact response of floor 81 contents. 
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(a) Initial impact configuration 
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(b) Calculated response of contents 
Figure 9-91. Calculated more severe WTC 2 impact response of floor 82 contents. 
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(a) Initial impact configuration 




(b) Calculated response of contents 
Figure 9-92. Calculated more severe WTC 2 impact response of floor 83 contents. 

9.7.2 Less Severe WTC 2 Global Impact Analyses 

Exterior Wall Damage for the Less Severe WTC 2 Impact 

A comparison of the south exterior wall observed and calculated damage from the less severe WTC 2 
global impact analysis is shown in Figure 9-93. The calculated impact damage to the exterior wall is 
shown with color fringes representing plastic strain magnitude, with undamaged sections in blue and 
strains at or above 5 percent shown in red. The wall region shown in Figure 9-93 includes column 402 to 
column 446, extending from floor 76 to floor 86 (spandrels at floors 77 through 86). The mode and 
magnitude of the calculated and observed impact damage on the exterior wall are still in good agreement 
in this less severe impact analysis. 
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(b) Calculated damage (t=0.66 s) 
Figure 9-93. Less severe impact damage to the WTC 2 exterior wall. 
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Core Structural Damage for the Less Severe WTC 2 Impact 

The overall model of the WTC 2 core structure and calculated response for the less severe global impact 
analysis is shown in Figure 9-94. The figure shows that the core had significant damage in the region 
close to the impact point. The columns in line with the aircraft fuselage failed on the impact side, and 
several of the core beams were also severely damaged or failed in the impact zone. In some cases, failure 
of the column splices located on floors 77, 80, and 83 contributed significantly to the failure of the core 
columns. 

The calculated damage to the core columns by row is shown in Figure 9-95, and the floor-by-floor 
damage to the core framing is shown in Figure 9-96. The core structures are shown with color fringes 
representing plastic strain magnitude, with undamaged sections in blue and strains at or above 5 percent 
shown in red. A summary of the column damage is provided in Table 9-12 and shown graphically in 
Figure 9-97, with the qualitative classification of the column damage levels provided previously in 
Figure 9-12. A total of three columns were severed, and two columns heavily damaged, compared to five 
severed columns and four heavily damaged columns in the base case WTC 2 impact analysis. 

Floor Truss and Slab Damage for the Less Severe WTC 2 Impact 

An overall front view of the floor truss structure in the WTC 2 impact zone, along with the calculated less 
severe impact damage to the trusses, is shown in Figure 9-98. The truss floor system on floors 79 
through 82 had sufficient damage from the impact that portions of the truss floor sections sagged 
downward as a result of the impact. A plan view of the calculated damage to the trusses on each floor is 
shown in Figure 9-99. The figure shows that the trusses experienced significant damage in the impact 
zone, with the heaviest damage on floor 81. The calculated impact response produced severe damage to 
the truss structures in the primary path of the fuselage. The truss structures were completely destroyed 
along the impact path on floor 81 from the exterior wall to the core. 
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(a) Initial geometry 




(a) Calculated impact damage 
Figure 9-94. Less severe impact response of the WTC 2 core. 
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Figure 9-95. Less severe impact response of the WTC 2 core columns. 
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(a) Floor 78 core framing damage 
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(b) Floor 79 core framing damage 




(d) Floor 81 core framing damage 




(e) Floor 82 core framing damage (f) Floor 83 core framing damage 

Figure 9-96. Less severe impact damage to the WTC 2 core floor framing (plan view). 
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Table 9-12. Core column response for the less severe WTC 2 impact. 



Column 


Location 


Damage Level 


Lateral Deflection of 
Column Centerline (in.) 


Column 801 


Floor 79 


Light 




Column 901 


Floors 79-80 


Severed 




Column 902 


Floors 78-79 


Heavy 


24 


Column 903 


Floors 77-81 


Severed 




Column 1001 


Floors 77-83 


Severed 




Column 1002 


Floors 80-81 


Moderate 




Column 1003 


Floors 79-80 


Heavy 


16 


Column 1004 


Floor 80 


Moderate 
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Figure 9-97. Summary of core column damage for the less severe WTC 2 impact. 
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(b) Calculated damage (t=0.661 s) 
Figure 9-98. Less severe impact damage to the WTC 2 floor truss (front view). 

The magnitude of truss floor damage was very similar for the base case and less severe global impact 
analyses. The parameters used in the less severe global impact analysis would primarily contribute to a 
reduced damage magnitude for the tower structures. However, the downward impact trajectory angle was 
increased from the 6 degree angle in the base case analysis to an 8 degree angle in the less severe impact 
analysis. This would have the effect of directing more of the impact energy downward, increasing the 
normal force on the floor structures in the impact zone. As a result, the combined effects of the analysis 
parameter variations produced very similar damage to the truss structure. 

The calculated damage to the WTC 2 floor slab for floors 78 through 83 for the less severe impact is 
shown in Figure 9-100. Similar to the damage levels in the floor truss structures, the magnitude of floor 
slab damage was very similar for the base case and less severe global impact analyses. The increase in 
the impact trajectory angle from 6 to 8 degrees offset the effects of the other analysis parameters for 
producing damage to the floor slab. 
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(c) Floor 81 truss damage 
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(d) Floor 82 truss damage 
Figure 9-99. Less severe impact damage to the WTC 2 floor truss (plan view). 
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(a) Floor 78 slab damage 




(c) Floor 80 slab damage 
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(e) Floor 82 slab damage 
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(f) Floor 83 slab damage 



Figure 9-100. Less severe impact damage to the WTC 2 floor slab (plan view). 
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Damage to WTC 2 Contents for the Less Severe Impact 

The calculated damage to the WTC 2 contents for the less severe impact is shown in plan views in 
Figure 9-101 through Figure 9-106 for floors 78 through 83, respectively. A comparison to the 
calculated damage for the base case WTC 2 impact analysis, shown previously in Figure 9-41 through 
Figure 9-46, indicated that the content damage zone was very similar, with a slight reduction in damage 
for the less severe impact. 
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(b) Calculated response of contents 
Figure 9-101. Calculated less severe WTC 2 impact response of floor 78 contents. 
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(a) Initial impact configuration 
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(b) Calculated response of contents 
Figure 9-102. Calculated less severe WTC 2 impact response of floor 79 contents. 
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(a) Initial impact configuration 
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(b) Calculated response of contents 
Figure 9-103. Calculated less severe WTC 2 impact response of floor 80 contents. 



324 



NISTNCSTAR 1-2B, WTC Investigation 



Global Impact Analyses 










II 








L i 


f i / I" 



mnmn 

LULJLUU 

■"UTLH 

LULJLUU 

rnnmnj^lnrnnrnnf^ 

LU LJ LU L [3 U LLTLJILTU [ 




(a) Initial impact configuration 




(b) Calculated response of contents 
Figure 9-104. Calculated less severe WTC 2 impact response of floor 81 contents. 



NISTNCSTAR 1-2B, WTC Investigation 



325 



Chapter 9 




nn 



LU U LU LJ |gl \l HI 






(a) Initial impact configuration 
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(b) Calculated response of contents 
Figure 9-105. Calculated less severe WTC 2 impact response of floor 82 contents. 
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(a) Initial impact configuration 



j j a 



I P n i 




(b) Calculated response of contents 
Figure 9-106. Calculated less severe WTC 2 impact response of floor 83 contents. 
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9.8.1 



COMPARISON OF DIFFERENT SEVERITY WTC 1 IMPACT ANALYSES 
Exterior Wall Damage for WTC 1 



The calculated damage to the north exterior wall of WTC 1 for the three different severity impacts is 
shown in Figure 9-107, along with a schematic of the observed damage. The overall agreement with the 
observed damage was good for all three analyses, with the base case global impact analysis providing the 
best match to the observed damage. The calculated damage magnitude was similar in each of the global 
analyses with small differences, as shown in the figure. 

The differences in apparent damage were largely due to panels that may have severed columns in one 
case and were removed at the connections in another. Visually this looks substantial as portions of panels 
that were damaged but still attached to the wall in the more severe impact, shown in the Figure 9-107(b), 
were completely removed in the less severe impact, shown in Figure 9-107(c). In general, the trend was 
for a larger opening produced by the less severe impact and a smaller opening in the more severe impact. 
The increase in the opening with reduced severity impacts can be explained by the increased tower panel 
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material toughness, resulting in a transition from column severing to connection bolt failure and panel 
removal. 
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(a) Base case impact damage 



(b) More severe impact damage 




(c) Schematic of observed damage 



(d) Less severe impact damage 



Figure 9-107. Impact damage to the north exterior wall of WTC 1. 

Toward the wing tips, where the columns and spandrels were not completely severed, the more severe 
impact damage analysis calculated higher damage to the exterior wall panels. These columns had the 
largest amount of material with plastic strains above 5 percent (shown in red in the figure). As would be 
expected, the less severe impact damage analysis calculated lower damage to the exterior wall, and the 
base case analysis calculated an intermediate level of damage near the wing tips. 



9.8.2 



Core Column Damage for WTC 1 



The calculated damage to the WTC 1 core columns for the three different severity impacts is shown in 
Figure 9-108. The figure shows that the core damage was concentrated in the region of the core closest 
to the impact point, and there was a clear correlation in damage magnitude with the impact severity. A 
total of one column was severed, and two columns were heavily damaged in the less severe impact, 
compared to three columns severed and four columns heavily damaged in the base case impact analysis 
and six columns severed and three columns heavily damaged in the more severe WTC 1 impact analysis. 
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(a) Calculated less severe impact analysis 




•Severed •Heavy OModerate 

(b) Calculated base case impact damage 
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(c) Calculated more severe impact analysis 
Figure 9-108. Comparison of core column damage for the WTC 1 impact analyses. 
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The strong correlation in the core damage with impact severity was expected. All of the parameter 
variations in the three increasing severity analyses would be expected to produce an increase in core 
damage. The flight parameters had an increasing impact speed and a shallower impact angle, directing 
more energy toward the core. The aircraft had an increasing weight and higher material toughness. The 
tower had reduced mass in the contents and a reduced material toughness. All of these variations 
contributed toward the increased core damage with increased impact severity. 

9.8.3 Floor Truss Damage for WTC 1 

The calculated damage to the WTC 1 floor truss structures for the three different severity impacts is 
shown in Figure 9-109. The magnitude of truss floor damage was very similar for the different severity 
WTC 1 global impact analyses. The parameters used in the increasing severity global impact analyses 
would primarily contribute to an increased damage magnitude for the tower structures. However, the 
downward impact trajectory angle was correspondingly reduced from 13.6 degrees in the less severe 
impact to 10.6 and 7.6 degrees in the base case and more severe impact scenarios, respectively. This 
would have the effect of directing less of the impact energy downward, reducing the normal force on the 
floor structures in the impact zone. As a result, the combined effects of the analysis parameter variations 
produced very similar damage to the truss structure for the three different severity impact scenarios. 
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(a) Calculated less severe impact analysis 
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(b) Calculated base case impact damage 



Column 
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(c) Calculated more severe impact analysis 
Figure 9-109. Comparison of floor truss damage for the WTC 1 impact analyses. 

9.9 COMPARISON OF DIFFERENT SEVERITY WTC 2 IMPACT ANALYSES 

9.9.1 Exterior Wall Damage for WTC 2 

The predicted impact damage to the south exterior wall of WTC 2 for the three different severity impacts 
is shown in Figure 9-110, with a schematic of the observed damage for comparison. As was the case for 
WTC 1, there were small differences in the predictions from each of the global analyses shown in the 
figure. Overall, the agreement with the observed damage from photographs was very good. The most 
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obvious differences were largely due to portions of panels that may have severed columns in one case or 
have been removed at the connections in another. 
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(a) Base case impact damage 
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(b) More severe impact damage 
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(c) Schematic of observed damage (d) Less severe impact damage 

Figure 9-110. Impact damage to the south exterior wall of WTC 2. 



9.9.2 



Core Column Damage for WTC 2 



The calculated damage to the WTC 2 core columns for the three different severity impacts is shown in 
Figure 9-111. The figure shows that the core damage was concentrated in the region of the core closest 
to the impact point, and there was a clear correlation in damage magnitude with the impact severity. A 
total of three columns were severed and two columns were heavily damaged in the less severe impact, 
compared to five columns severed and four columns heavily damaged in the base case analysis and ten 
columns severed and one column heavily damaged in the more severe WTC 2 impact analysis. 

The strong correlation in the core damage with impact severity was expected. All of the parameter 
variations in the three increasing severity analyses would be expected to produce an increase in core 
damage. The flight parameters had an increasing impact speed and a shallower impact angle, directing 
more energy toward the core. The aircraft had an increasing weight and higher material toughness. The 
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tower had reduced mass in the contents and a reduced material toughness. All of these variations 
contributed toward the increased core damage with increased impact severity. 
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(a) Calculated less severe impact analysis 
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(b) Calculated base case impact damage 
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(c) Calculated more severe impact analysis 
Figure 9-111. Comparison of core column damage for the WTC 2 impact analyses. 
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9.9.3 



Floor Truss Damage for WTC 2 



The calculated damage to the WTC 2 floor truss structures for the three different severity impacts is 
shown in Figure 9-112. The magnitude of truss floor damage was very similar for the different severity 
WTC 2 global impact analyses. The parameters used in the increasing severity global impact analyses 
would primarily contribute to an increased damage magnitude for the tower structures. However, the 
downward impact trajectory angle was correspondingly reduced from 8 degrees in the less severe impact 
to 6 and 5 degrees in the base case and more severe impact scenarios, respectively. This would have the 
effect of directing less of the impact energy downward, reducing the normal force on the floor structures 
in the impact zone. As a result, the combined effects of the analysis parameter variations produced very 
similar damage to the truss structure for the three different severity impact scenarios. 
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Floor 78 



(a) Calculated less severe impact analysis 
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(b) Calculated base case impact damage 
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(c) Calculated more severe impact analysis 
Figure 9-112. Comparison of floor truss damage for the WTC 2 impact analyses. 
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9.10 COMPARISON BETWEEN WTC 1 AND WTC 2 

The comparison of the aircraft impact response and resulting tower damage for WTC 1 and WTC 2 was 
complicated by the differences in the two impact scenarios. The base case WTC 1 impact was close to 
centered and perpendicular on the face of the tower, with the long span trusses between the impact point 
and the core. The base case WTC 1 impact scenario resulted in a debris trajectory where almost all of the 
aircraft debris passed through the core. The baseline impact conditions for WTC 1 were a 443 mph 
collision with a downward impact trajectory angle of 10.6 degrees. In contrast, the baseline WTC 2 
impact was off center and angled away from the core, resulting in a significant fraction of the aircraft 
debris cloud outside (east) of the core. The WTC 2 impact had short span trusses between the impact 
point and the core. Finally, the baseline impact conditions for WTC 2 were a 546 mph collision with a 
downward impact trajectory angle of 6 degrees. 

Exterior Wall Damage 

The calculated exterior wall damage for the base case WTC 1 and WTC 2 impacts are compared in 
Figure 9-113. Despite the differences in impact conditions, the mode and magnitude of damage to the 
exterior walls were quite similar in both towers. This was because the impact loads distributed over the 
majority of the aircraft structures were much larger than the exterior column rupture strength. The details 
of the failure mode (column deformation and rupture or failure and separation of bolted column end 
connections) were determined by the proximity of the floor slab and column joints to the impact point. 
For both impacts, the wing tip structures imparted damage, but did not completely fail the columns. 

Core Column Damage 

The calculated core column damage for the base case WTC 1 and WTC 2 impacts are compared in 
Figure 9-114. In the WTC 1 impact, there were three columns severed and four columns heavily 
damaged. The calculated region of significant core column damage appears to extend three column rows 
deep into the core. In contrast, the calculated damage for the WTC 2 impact included five columns 
severed and four columns heavily damaged, and the region of significant core column damage appears to 
extend four column rows deep. This increase in core damage was even more significant since the impact 
zone was 15 floors lower in WTC 2 (and therefore designed to carry more gravity loads), and as a result 
the core columns were heavier and more resistant to impact damage in the WTC 2 impact zone. 

The differences in the core column damage between WTC 1 and WTC 2 can be explained by two primary 
factors. The first was that the WTC 2 impact speed was 23 percent higher (approximately 50 percent 
higher impact energy), and the shallower impact angle directed more impact energy inward toward the 
core. The second factor was that the orientation of the core relative to the impact was different in the two 
towers, as the core was closer to the impact point in WTC 2. As a result, WTC 2 had reduced energy 
absorbing capacity due to the shorter floor structures and less building contents between the impact point 
and the core. 
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(a) WTC 1 calculated damage (t=0.715 s) 
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(b) WTC 2 calculated damage (t=0.62 s) 
Figure 9-113. Comparison of base case impact damage to the exterior wall. 
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(a) WTC 1 calculated damage (t=0.715 s) 
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(b) WTC 2 calculated damage (t=0.62 s) 
Figure 9-114. Comparison of base case impact damage to the core columns. 

Floor Truss Damage 

The calculated floor truss damage for the base case WTC 1 and WTC 2 impacts are compared in 
Figure 9-115. The comparison shows that the WTC 1 floor truss had greater damage and collapse of the 
truss floor despite the lower aircraft impact energy. The higher amount of truss floor damage and 
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deflection in WTC 1 can be explained by two factors. The primary factor was that the WTC 1 downward 
impact trajectory was nearly twice as steep as that of the WTC 2 impact. As a result, the steeper impact 
angle directed more impact energy normal to the floor slab. The vertical component of the impact load in 
WTC 1 was approximately 40 percent higher than in WTC 2 (the vertical momentum of the aircraft that 
impacted WTC 1 was about 40 percent higher than that for the aircraft that impacted WTC 2). The 
secondary factor was that the damage to the long span truss floors in the WTC 1 impact zone produced 
larger displacements than a corresponding damage level to a short span truss region. 
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(a) WTC 1 calculated damage (t=0.715 s) 
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(b) WTC 2 calculated damage (t=0.715 s) 
Figure 9-115. Comparison of base case impact damage to floor truss (front view). 



9.11 



COMPARISON WITH OBSERVABLES 



The observable evidence available to help validate the global impact analyses included the following: 

• Damage to the buildings' exteriors documented by photographic evidence. 

• Floor damage visible from the buildings' exteriors documented by photographic evidence. 

• Aircraft debris external to the towers as documented by photographic evidence. 

• Eyewitness accounts from survivors who were inside portions of the buildings. 
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The most valuable observable from a modeling standpoint was the damage to the impacted exterior wall 
of each tower. The impact damage to the exterior tower walls was well documented, and the response did 
not depend much on unknown parameters, such as the detailed office layout on each floor. Good 
agreement of the calculated and observed damage profile indicated that the geometric modeling of the 
aircraft and the initial trajectory and orientation of the aircraft were accurate. The agreement of both the 
mode and magnitude of the structural damage on the impact wall served to partially validate the 
constitutive and damage modeling of the aircraft and exterior wall structures of the tower. The agreement 
in exterior wall damage, based on the modeling methodologies described in this report, contributed to the 
confidence that the damage predictions for the interior of the towers were reasonably estimated. 

The calculated and observed impact damage to the exterior walls were described previously in this 
chapter and are summarized in Figure 9-116 for the base cases. There were small differences in the 
damage details, but the overall agreement with the observed damage is very good. 
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(a) Base case WTC 1 impact analysis 



(b) Schematic of observed WTC 1 damage 




(c) Base case WTC 2 impact analysis (d) Schematic of observed WTC 2 damage 

Figure 9-116. Impact damage to the exterior walls of the WTC towers. 
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9.11.1 Comparison with Observables on WTC 1 

Damage on the South Exterior Wall of WTC 1 

The exterior panel from column 329 to column 331 between floors 94 through 96 on the south face of 
WTC 1 was knocked free by landing gear and possibly other debris. These columns were located in the 
center of the south wall of WTC 1, as shown in Figure 9-117. In both the base case and more severe 
damage global analyses, aircraft debris impacted the south face of the tower, as shown in Figure 9-118 
and Figure 9-119, and exited the building. In the less severe damage analysis, as shown in Figure 9-120, 
none of the aircraft debris that passed through the core was calculated to exit the building. The figures 
also show the calculated landing gear debris for all simulations. None of the debris impacting the south 
wall happened to contain landing gear fragments. In the base case analysis, the debris impacted columns 
328 to 330 at floor 96. In the more severe impact analysis, debris impacted columns 328 to 333 on both 
floors 95 and 96. In the base case analysis, very little damage was done to the exterior panels on the 
south wall. However, damage was heavy in the more severe damage analysis, as shown in Figure 9-121. 
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Figure 9-117. Representative floor plan (based on floors 94 to 95 of WTC 1). 
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(a) Calculated aircraft debris (t = 0.71 s) 
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(b) Calculated landing gear debris (t = 0.71s) 
Figure 9-118. Base case aircraft debris distribution in WTC 1. 
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(a) Calculated aircraft debris (t = 0.685 s) 
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(b) Calculated landing gear debris (t = 0.685 s) 
Figure 9-119. More severe damage aircraft debris distribution in WTC 1. 
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(b) Calculated landing gear debris (t = 0.735 s) 
Figure 9-120. Less severe damage aircraft debris distribution in WTC 1. 
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Figure 9-121. Damage to the south face of WTC 1 from the more severe damage 

global analysis. 

Because of model size constraints, the panels on the south face of WTC 1 were modeled with a very 
coarse resolution. Neither the spandrel splice joints nor exterior column butt joints were modeled. 
Column ends and spandrel edges were merged together. The model, therefore, underestimates the 
damage to the tower on this face. The calculated damage produced by the more severe impact is shown in 
Figure 9-121. Columns 329-331 on floors 94 through 96 had sustained substantial damage. Had a fine 
mesh been used on these columns, it is likely that they would have failed on floor 95, and possibly on 94 
and 96. Based on the failure modes observed on the north face and on the speed and mass of the debris, 
the panel would potentially be knocked free by failing at the connections. 

Landing Gear Trajectory 

It is believed that a portion of the main landing gear of AA 11 exited WTC 1 at the 94th or 95th floor and 
landed at the corner of Rector and West Streets. This debris is believed to be a tire, wheel, brake 
assembly, and hub of a main landing gear, as shown in Figure 9-122. Based on the final position of the 
landing gear and assuming the landing gear to be a projectile with a horizontal initial velocity, the exit 
speed of the landing gear from the south wall of WTC 1 can be estimated to be about 105 mph. Note that 
there is a significant uncertainty in this estimate associated with the exit trajectory, aerodynamic effects, 
landing position rather than final resting position of debris, etc. Another piece of landing gear debris, 
shown in Figure 9-123, was found embedded in what is postulated to be the panel containing 
columns 329, 330, 331, running from the 93rd to the 96th floors. This panel was dislodged from the 
building and found at Cedar Street near its intersection with West Street. As little other damage had been 
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documented on the south face of WTC 1, it is postulated that the landing gear debris that landed at the 
corner of Rector St. and West St. also exited through this panel location. 

The amount of aircraft debris found to exit WTC 1 in the global impact analyses varied, as shown in 
Figure 9-118 to Figure 9-120. However, no portion of the landing gear was observed to exit the tower in 
the simulations, but rather was stopped inside, or just outside, of the core. In order to simulate the 
trajectory of specific pieces of aircraft debris, a fairly precise knowledge of the internal configuration of 
the building would be needed. This is especially true with components passing through the core of the 
building, where some of the most massive building contents and partition walls were present. 
Uncertainties regarding the internal layout of each floor, such as the location of hallways or walls, can 
make the difference between debris from a specific component passing through or being stopped inside 
the structure. Modeling uncertainties may also have contributed to the inability to predict the trajectory of 
specific aircraft components. 




Figure 9-122. Landing gear found at the corner of West and Rector streets. 
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\© 2001 New York City Police Department. All rights reserved! 
Figure 9-123. Landing gear found embedded in exterior panel knocked free from WTC 1. 



Stairwell Disruption 

According to eyewitness interviews, stairwells 1 (referred to also as stairwell A), 2 (stairwell C), and 3 
(stairwell B) inside the core were impassable at floor 92 and possibly above after the impact of AA 11. 
The calculated base case stairwell disruption is shown in Figure 9-124 for floors 93 through 97. Stairwell 
positions are outlined with red boxes in the figure. No debris or disruption was observed to the core on 
floor 92 in the calculation; therefore, it is not shown in the figure. Recall that the global model for 
WTC 1 only contained partition walls in the core on floors 94 through 97. Therefore, the ability to 
ascertain damage and/or debris in the stairwell on floors 92 and 93 was limited. The floor slab was 
removed from view on floors 94 through 97 so that debris is more visible. 

Based on the calculated damage to, or debris in, the stairwells on floors 94 to 96, all three stairwells 
appear impassable. Given that falling debris in these areas would cause further subsequent damage to the 
floors below, as well as block passage on these floors, this result was reasonably consistent with the 
eyewitness accounts. 

Floor Damage Visible on the North Face of WTC 1 

One location where the damage to the WTC 1 truss floors could be observed was through the opening in 
the tower exterior produced by the aircraft impact. A photograph of the impact damage on the north face 
of WTC 1 is shown in Figure 9-125(a). The magnitude of damage is difficult to quantify as a result of 
the strong contrast in lighting between the tower interior and exterior and the smoke inside the building. 
However, the photograph shows that the truss floor was heavily damaged and/or removed in the primary 
impact zone. The depth of the floor damage extending into the tower could not be determined. 
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(d) Floor 96 



(e) Floor 97 



Figure 9-124. Base case stairwell disruption in WTC 1. 
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(a) Damage photograph 




(b) Calculated damage 
Figure 9-125. Calculated and observed WTC 1 damage (front view). 
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A corresponding image of the calculated damage to the tower structures is shown in Figure9-125(b). The 
structures beyond the start of the core were removed and replaced with a black background for 
comparison with the photograph. Although a quantitative comparison of the calculated and observed 
damage could not be made from the available damage photographs, the truss floor damage appeared to be 
consistent. 

9.11.2 Comparison with Observables on WTC 2 

Damage on the North Exterior Wall 

From photographic evidence, such as that shown in Figure 9-126, damage on the north wall at the 
northeast corner of WTC 2 was documented and is shown in Figure 9-127. As mentioned earlier, there 
was significant uncertainty as to the actual layout of the workstations and other building contents on the 
impacted floors of the towers. Recall that generic workstation configurations were used to model these 
building contents, as shown in the northeast corner of WTC 2 in Figure 9-128(a). Uncertainties regarding 
this layout, such as missing partition walls and workstations, can make the difference between debris 
from a specific component passing through or being stopped inside the structure. The base case impact 
response of the northeast corner of WTC 2 on the 81st floor is shown in Figure 9-128(b). 




Figure 9-126. Impact Damage to the Northeast Corner of the Exterior Wall of WTC 2. 
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Figure 9-127. Documented damage to the 81st floor of the northeast corner of WTC 2. 
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(c) Residual speed of aircraft debris (contours in mph) 
Figure 9-128. Base case analysis on the 81st floor of the northeast corner of WTC 2. 
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Aircraft debris on the 81st floor is shown in Figure 9-128(c), with the coloring depicting the residual 
speed of the debris field. Notice that some of the debris in this figure, weighing approximately 3,800 lb, 
was traveling at 110-150 mph and was projected to impact between columns 252 and 256. The leading 
debris was in fact portions of the starboard main landing gear main strut and main landing gear beam. 
That significant debris was projected to impact in the region of significant damage shows positive 
agreement with damage evidence available for the north wall of WTC 2. 

Stairwell Disruption 

According to eyewitness interviews, stairwells 2 and 3 on floor 78 of WTC 2 were impassable. 
Stairwell 1 (referred to also as stairwell A), which was located in the northwest corner of the core, was 
passable. The calculated base case stairwell disruption is shown in Figure 9-129. Stairwells 1 and 2 
(stairwell C) on floor 78 of WTC 2 were outside of the core column region. These stairwells were not 
included in the WTC 2 model. Therefore, a good assessment could not be made for stairwell 2. 
However, disruption to stairwell 3 (stairwell B) is shown in Figure 9-129. From the damage shown in the 
figure, the stairwell appears to be impassable. As no damage or debris was seen in the northwest corner 
of the core, the top right in the figure, stairwell 1 in this area of the core was likely unaffected. Both of 
these assessments were consistent with the eyewitness accounts. 
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Figure 9-129. Base case stairwell disruption on floor 78 in WTC 2. 

Landing Gear Trajectory 

A portion of the landing gear of UAL 175 exited WTC 2 and landed on the roof of 45 Park Place. No 
photographic evidence was available to document the size of the fragment and whether this was a nose or 
main landing gear. From the damage to the building, it was believed that the landing gear fragment might 
have exited somewhere along the north wall between column 251 and the northeast corner on floor 81. 
Based on the final position of the landing gear and assuming the landing gear to be a projectile with a 
horizontal initial velocity, the exit speed of the landing gear from the north wall of WTC 2 can be 
estimated to be about 102 mph. Note that there is a significant uncertainty in this estimate associated with 
the exit trajectory, aerodynamic effects, landing position rather than final resting position of debris, etc. 
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The calculated aircraft debris distribution and landing gear and engine debris distributions for UAL 175 
are shown in Figure 9-130 to Figure 9-132 for the three severity scenarios. A portion of the port main 
landing gear was seen to exit the building at approximately 230 mph in the more severe impact analysis, 
as shown in Figure 9-131(b). No landing gear debris exited the building in either the base case or less 
severe simulations. At the conclusion of the simulation, the base case analysis had a substantial piece of 
the starboard main landing gear still at approximately 130 mph that was expected to impact the northeast 
corner. 

Engine Trajectory 

A portion of an engine also exited the tower at the northeast corner of the building and was found at the 
intersection of Murray and Church Streets. From the damage to the building, it was believed that the 
engine exited the building in this corner of WTC 2. Based on this trajectory, it was estimated that the 
engine exited the building at approximately 120 mph. The engine trajectories predicted from the base 
case global analysis are shown in Figure 9-133, which indicates that the engine that exited from the 
northeast corner is likely the starboard engine. The dotted line indicates the extrapolated engine flight 
path based on the initial trajectory of the starboard engine. Notice that this trajectory would result in 
engine fragments exiting at the northeast corner. In the simulations, the engines were projected to stop 
short of this position, although they followed the extrapolated trajectory reasonably well. 

Speed histories for the aft portion of the starboard engine are shown in Figure 9-134. The engine would 
typically break up into smaller fragments from the forward section of the engine and a larger section from 
the aft end, as shown in Figure 9-135. In all three simulations, the speed was seen to drop by 
approximately 200 mph due to impact with the exterior panel, floor slab, and floor truss. Interaction with 
these portions of the structure ends by approximately 0.12 s. This initial impact is shown in 
Figure 9-136. The engine debris then continued through the tenant space of the 81st floor, plowing 
through the workstations and contents. Whether or not the fragment passes over these contents, or if 
other debris and fuel remove the contents from the engine's path, affects the deceleration of the fragment. 
At the end of the simulation, the speed of the aft portion of the engine was below 80 mph, and it was more 
than 60 ft from the northeast corner of the building. For these calculations, it was estimated that the 
building contents would likely stop the engine fragment prior to impacting the northeast corner of the 
exterior wall. 

None of the three WTC 2 global impact simulations resulted in a large engine fragment exiting the tower. 
However, the impact behavior suggests that only minor modifications would be required to achieve this 
response. For example, if the starboard engine impact location was lowered by 1 to 2 ft, which is within 
the aircraft impact geometry uncertainty range, the engine would likely have a greater residual speed 
inside the tower (over 100 mph). In the global analyses performed, the engine impacted the underside of 
the 82nd floor, as shown in Figure 9-136. This resulted in a large reduction in speed of approximately 
200 mph. In the component analyses discussed in Chapter 5, the engine speed decreased by roughly 
60 mph when impacting an exterior panel alone. This additional speed would likely result in a large 
engine fragment exiting the northeast corner of the tower. 
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(b) Calculated engine and landing gear debris (t = 0.62 s) 
Figure 9-130. Base case damage aircraft debris distribution in WTC 2. 
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(b) Calculated engine and landing gear debris (t = 0.58 s) 
Figure 9-131. Aircraft debris distribution in the more severe WTC 2 impact. 



NISTNCSTAR 1-2B, WTC Investigation 



355 



Chapter 9 



It y a n ncmat a i p a ii 







•.•US'' 






'*.;** 






$ 



^ * 



L, 



■gpn- a iaigpg B dPfl p cm 



(a) Calculated aircraft debris (t = 0.66 s) 



■j a a | g a y a no il 8 ii t o 



a B a B □ n -g g □ g y j 6 n-a- g a □ a a a B B 6 n 



»V 



;*+.., 



*..'* 



> ..,-■ ' 



(b) Calculated engine and landing gear debris (t = 0.66 s) 
Figure 9-132. Aircraft debris distribution in the less severe WTC 2 impact. 
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Figure 9-133. Starboard engine fragment trajectory in the base case global analysis of 

WTC 2. 
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Figure 9-134. Speed of the aft portion of the starboard engine. 
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(a) Undamaged engine 

2001 GeoTge Marengo 



(b) Large engine fragment 




(c) Engine fragment found at Murray and Church St. 
Figure 9-135. Calculated and observed engine damage. 
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Figure 9-136. Starboard engine impact with the south face of WTC 2 in the base case 

global analysis. 

Other minor modifications to the model could also result in a large engine fragment exiting the building. 
As mentioned previously, there was significant uncertainty in the distribution of building contents on the 
floors of the impact area. If any portion of the east side of WTC 2 was relatively free of office materials, 
the engine fragment would have been free to move relatively unrestricted and would have experienced 
little loss of speed. After the engine entered the structure, and without office materials, the engine 
fragment would only slow due to friction with the floor slab and occasional interaction with floor trusses 
above. After initially entering the building, the engine did not further penetrate the floor slab. Removing 
much of these building contents from the east side would result in the starboard engine fragment 
impacting the northeast corner of the tower with sufficient speed to exit the building. Little or no 
difference in core damage would result, as debris in this area had no chance of impacting the core. 

Floor Damage Visible on the South Face of WTC 2 

One location where the damage to the WTC 2 truss floor could be observed was through the opening in 
the tower exterior produced by the aircraft impact. A similar comparison for the WTC 1 truss floor 
damage was shown in Figure 9-125. The magnitude of damage was difficult to quantify as a result of the 
strong contrast in lighting between the tower interior and exterior and the smoke inside the building. This 
was worse for WTC 2, where the prevailing wind and fire conditions resulted in larger quantities of 
smoke exiting through the opening on the impact face. The partial photographic evidence did suggest that 
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a similar level of truss floor damage in the impact zone occurred for WTC 2. The severity and the depth 
of the floor damage extending into the tower could not be determined. Although a quantitative 
comparison of the calculated and observed damage could not be made from the available damage 
photographs, the truss floor damage appears to be consistent. 

The 'Cold Spot' on the North Face of WTC 2 

A 'cold spot' was observed on the north face of the tower between columns 238 and 250 on floors 80, 81, 
and 82. The cold spot was a region of the tower where no debris could be seen from the exterior of the 
tower, and no significant fires were observed prior to tower collapse. 

Much of the explanation for the cold spot was obtained from an analysis of the debris trajectory aligned 
with the cold spot. The debris path, obtained by projecting the width of the cold spot along the initial 
lateral impact trajectory of the aircraft, is shown in Figure 9-137 (13 degrees relative to the tower face 
normal). This region was aligned laterally with the left side of the fuselage and the port wing structures. 
Considering the baseline impact orientation and trajectory, shown in Figure 9-138, it can be seen that 
much of the wing debris impacted on floors lower than the observed cold spot. Only debris from very 
close to the fuselage would be expected on floor 80 or above. The debris from the port wing, including 
the majority of the aircraft fuel in the left side tanks, entered at floors 78 and 79. 

The base case WTC 2 global analysis calculated a small amount of aircraft debris passing through the 
cold zone on floors 80 and 81. However, the building contents were not completely modeled over the 
entire path in this section. After clearing the core region, the debris in the calculation had primarily an 
open path to the cold spot on the North wall of WTC 2. If all of the internal contents had been included, 
it is likely that all of this debris would have been stopped before reaching the cold spot. 

The comparison of the calculated and observed impact response cold spot is inconclusive. Much of the 
absence of damage and aircraft debris in this region is explained by the impact orientation and trajectory. 
Much of this region is not directly in the path of significant aircraft fuel and debris. In addition, the debris 
aligned with the cold spot would be required to pass through a significant portion of the core. A more 
accurate analysis of the impact mechanics leading to the formation of a cold spot would require a specific 
survey of the tenant layout, including both contents that acted as a barrier to the debris and walls that 
provided a barrier to subsequent fire propagation. 



NISTNCSTAR 1-2B, WTC Investigation 361 



Chapter 9 



/~~~~ 



|l I I « - | — r- p | a f ~T~~ i < | I H 

!= = = = l[--i-j=i-i II I- j= = = = H 

*- ~ - " - — i-1 ~ ~ ~i I ' 

i = = '^^^^^ = - f ■ 

« = = =*c^==«-^ i V ^i * I rp = = = = =. 
n I r I T "tFw ~3 n V I I I *■ I ' P 



Ti nOnODODODaDCDCDCnCIDODODQDCDCDOnDBCDQDQDQDflDODCnQDODODQD l D r 

Figure 9-137. Projected debris path for the WTC 2 north face cold spot. 
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Figure 9-138. Base case WTC 2 impact orientation and trajectory 
(vertical approach angle = 6°, lateral approach angle = 13°). 
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9.11.3 Comparison with Observables - Summary 

A summary of the observables, their relative significance, and the level of agreement with these 
observables seen in the simulations is shown in Table 9-13 and Table 9-14 for WTC 1 and WTC 2, 
respectively. The most valuable observable was the damage to the impacted exterior wall of each tower. 
The impact damage to the exterior walls was well documented and the response did not depend as much 
on unknown parameters, such as the detailed office layout on each floor. Eyewitness accounts of damage 
to stairwells and visible floor damage were also significant as they were the only available data on the 
damage to the tower interiors. That these observables were in good agreement adds greater credibility to 
the predicted damage to the tower interiors. 

Observed trajectories of specific aircraft components, such as the landing gear and engines, were 
considered to be of lower importance in validating the simulated damage to the tower. A fairly precise 
knowledge of the internal configuration of the building would be needed in order to simulate the 
trajectory of specific aircraft debris. Damage to the opposite side of each tower from the point of impact 
was also of lower importance. These parts of the tower were modeled with lower resolution and as a 
result, the models were not sufficient to capture the detailed damage. The actual damage was also 
predominantly caused by the landing gear and engine components already discussed. Overall, agreement 
between the observables and simulations is good. 



Table 9-13. Comparison with observables from WTC 1. 



Observable 


Importance 


Agreement 


Damage to the north exterior wall 


Very significant 


Very good 


Damage to the south exterior wall 


Slightly significant 


Fair 


Landing gear trajectory 


Slightly significant 


Poor 


Stairwell disruption 


Significant 


Good 


Floor damage visible on the north 
face 


Significant 


Good (unclear) 



Table 9-14. Comparison with observables from WTC 2. 



Observable 


Importance 


Agreement 


Damage to the south exterior wall 


Very significant 


Very good 


Damage to the north exterior wall 


Slightly significant 


Good 


Stairwell disruption 


Significant 


Good 


Landing gear trajectory 


Not significant 


Fair 


Engine trajectory 


Slightly significant 


Good 


Floor damage visible on the south 
face 


Significant 


Good (unclear) 


The 'cold spot' on the north face 


Slightly significant 


Unclear 



9.12 



COMPARISON WITH PREVIOUS STUDIES 



Two previous studies were conducted to estimate the impact damage to the WTC towers. These studies 
were performed by staff from Massachusetts Institute of Technology (MIT) (Wierzbicki, Xue, and 
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Hendry-Brogan 2002) and Weidlinger Associates, Inc. (WAI) (Levy and Abboud 2002). The MIT study 
used an energy balance approach to estimate damage to the core columns. Estimates were made for the 
initial kinetic energy of the impacting aircraft and the internal energy absorbed in fragmentation of the 
aircraft and damage to the tower exterior columns, floor slab, and core columns. The energy absorbed by 
the core was used to estimate the number of failed core columns. 

The WAI study used the FLEX finite element code to calculate the aircraft impact damage to both towers. 
The FLEX family of finite element modeling software (Vaughan 1997) was developed and maintained by 
WAI. FLEX is an explicit, nonlinear, large deformation transient analysis finite element code for the 
analysis of structures subjected to blast, impact, and shock loadings. The overall code architecture is very 
similar to that of LS-DYNA, used to calculate the aircraft impact damage in this investigation. 

In the WAI calculations, the aircraft and WTC towers models were composed of beam and shell elements. 
The aircraft model consisted of 27,000 shell elements and 23,000 beam elements. The aircraft fuel was 
included in the model by increasing the mass of the structures in the wing box. The tower models 
included the exterior wall on the impact face, the floor structures, and the core frame for floors 91-101 
and floors 76 through 86 for WTC 1 and WTC 2, respectively. The tower models had fixed boundary 
conditions at the top and bottom floors. 

Comparison of Exterior Wall Damage 

The calculated base case impact damage to the exterior north wall of WTC 1 from this study is compared 
to the impact damage calculated by WAI in Figure 9-139. The figure also shows a schematic of the 
damage observed in photographic evidence. Figure 9-140 shows a similar comparison for the south wall 
of WTC 2. In both towers, the base case impact damage estimated in this study closely matched the 
observed damage. The damage profiles in the WAI impact simulations had some noticeable differences. 
The first is that the damage profile included complete failure of the exterior columns over the entire 
length of the wings and to the top of the vertical stabilizer. The second difference is that the failure mode 
of the exterior walls was dominated by local rupture of the columns adjacent to the impact point, with less 
influence of the bolted connections on panel failure and removal. 

The differences in the damage profiles in the two calculations most likely resulted from a variety of 
differences in the models. One major difference between the two studies was in the fidelity of the aircraft 
models. The WAI Boeing 767 model was based on their model of a Lockheed C-141B military transport. 
In the WAI model, the external geometry of the C-141B was modified to fit the dimensions of the 
Boeing 767, but the internal components, such as stiffener configuration and material thicknesses and 
properties, remained the same. The differences in the internal structure and materials could affect the way 
the aircraft responded to the impact. The aircraft model used in this study also contained an order of 
magnitude more elements (70,000 bricks, 562,000 shells, and 61,000 SPH particles) than the WAI model 
(27,000 shell elements and 23,000 beam elements). The higher resolution of the NIST model could also 
account for significant differences in the determination of the impact load distribution and resulting 
exterior damage. Additionally, the NIST model explicitly modeled the fuel. If the fuel mass in the WAI 
model was spread out further toward the wing tips as part of the wing structure, it would be expected that 
the calculated column damage would extend over a wider portion of the wings. 
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(a) Schematic of observed damage 
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(b) NIST base case impact damage 
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(c) WAI calculated damage (Levy and Abboud, 2002) 
Figure 9-139. Comparison of impact damage to the WTC 1 exterior wall. 
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(a) Schematic of observed damage 
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(b) NIST base case impact damage 




c) WAI calculated damage (Levy and Abboud, 2002) 
Figure 9-140. Comparison of impact damage to the WTC 2 exterior wall. 

Secondary differences in the WAI and NIST impact analyses included, but were not limited to, variations 
in impact conditions (impact velocity, angle, location, etc.), aircraft model differences (airframe 
geometry, component thicknesses, mass distribution, material properties, etc.) and tower model 
differences (material properties, geometry, joint modeling, number of elements, etc.). 
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Comparison of Core Column Damage 

Table 9-15 compares the estimated core column damage from the various studies. For WTC 1, MIT 
(Wierzbicki, Xue, and Hendry-Brogan 2002) estimated that four to 12 core columns failed. This MIT 
estimate of core columns was based on energy balance calculations and corresponded to a damage 
distribution ranging from four columns failed over a three-story length to 12 columns failed over a single 
floor length. The expected distribution of damage would fall between these bounds, with some columns 
damaged on a single floor and others with damage distributed on multiple floors. WAI gave two 
estimates for core column failure. The first estimate of 23 core columns failed and five damaged was 
obtained from the FLEX impact analysis. The second estimate of 20 failed columns was the number used 
in their collapse analysis. The NIST base case impact damage of three severed and four heavily damaged 
and less severe estimate of one severed and two heavily damaged falls below both the MIT and WAI 
estimates. The more severe estimate of six severed and three heavily damaged falls in the middle of the 
MIT range, but still well below the WAI estimates. 



Table 9-15. Comparison of damage to core columns from various studies. 



WTC Impact 
Investigation 


WTC 1 Core Column Damage 


WTC 2 Core Column Damage 


MIT 

Impact Analysis 


4-12 Failed 


7-20 Failed 


WAI 

Impact Analysis 


23 failed and significantly damaged 
Plus 5 damaged 


14 failed and significantly damaged 
Plus 10 damaged 


WAI 

Collapse Analysis 


20 Failed 


5 Failed 


NIST Base Case 
Impact Analysis 


3 Severed 
Plus 4 Heavily Damaged 


5 Severed 
Plus 4 Heavily Damaged 


NIST More Severe 
Impact Analysis 


6 Severed 
Plus 3 Heavily Damaged 


10 Severed 
Plus 1 Heavily Damaged 


NIST Less Severe 
Impact Analysis 


1 Severed 
Plus 2 Heavily Damaged 


3 Severed 
Plus 2 Heavily Damaged 



A similar trend in the predicted damage to the core columns was found in the WTC 2 analysis. MIT 
estimated seven to 20 columns failed (from seven columns failed over a three-story length to 20 columns 
failed over a single floor length). WAI calculated 14 core columns failed and another 10 damaged in their 
FLEX analysis, but reduced the number of failed columns to five for their collapse analysis. The NIST 
base case impact damage of five severed and four heavily damaged, as well as the more severe estimate 
of 10 severed and one heavily damaged fall in the middle of the range predicted by MIT. The less severe 
impact scenario predicted fewer columns severed and heavily damaged than the MIT and WAI studies. 

The MIT prediction for the number of failed core columns agrees remarkably well with the NIST 
estimates using their simplified analysis. Differences may be a result of the estimates of material 
properties and structural geometry used (MIT did not have access to the detailed structural drawing of the 
WTC towers for their study), approximations in the estimates of damage mode and resulting energy 
absorption, as well as the fact that the MIT study did not include the energy absorbed by internal tower 
contents. 
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The WAI impact analysis predicted much higher core column failure and damage than the NIST 
estimates. One reason for the higher damage prediction may be attributed to the lack of internal tower 
contents, such as workstations and other live loads, in the WAI model. This study found that the internal 
tower material absorbed a significant amount of the impact energy and, therefore, reduced the loads 
applied to the core columns. Another reason for the higher damage prediction in the WAI study could 
result from the aircraft model. As noted above, the WAI aircraft impact simulation overpredicted the 
extent of column damage and failure on the exterior wall. It is possible to assume that the aircraft model 
would also overpredict the damage to the core columns, especially that this damage configuration resulted 
in an unstable tower (Levy and Abboud, 2002). 

In conducting a collapse analysis, WAI used engineering estimates to reduce the number of failed 
columns from that predicted by their FLEX model to stabilize the tower immediately after impact. 
Despite this adjustment, the WAI study still estimates significantly higher damage for WTC 1 than the 
MIT and NIST studies. For WTC 2, their adjusted estimate falls in line with the MIT and NIST studies. 

In general, the MIT and WAI studies appear to over-predict the damage to the core columns compared to 
the NIST estimates. 
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10.1 INTRODUCTION 

Approximate, simplified impact analyses were carried out to provide guidance to the global finite element 
impact analyses. The primary purpose of these approximate analyses was to support the finite element 
impact analyses and identify aspects of the impact physics that needed to be included, or could be 
eliminated from the models. In this chapter, a few of the supporting approximate analyses are presented. 
The specific analyses described include: 

• The analysis of the overall aircraft impact forces and assessment of the relative importance of 
the airframe strength and weight distribution. 

• The evaluation of the potential effects of the energy in the rotating engine components on the 
calculated engine impact response. 

• The influence of the static preloads in the towers on the calculated impact damage and 
residual strength predictions. 

• The analysis of the load characteristics required to damage core columns compared to the 
potential loading from impact of aircraft components. 

10.2 SIMPLIFIED ANALYSIS OF AIRCRAFT IMPACT 

A simple analysis of the aircraft impact was performed at the initiation of this study to quantify the 
relative importance of the airframe structural strength versus the effects of the nonstructural mass and 
inertial effects on the impact loads and response. The primary objective of this analysis was to develop 
an improved understanding of the impact physics prior to development of the global impact models. 

Most prior aircraft impact analyses were concerned with the collision against nuclear reactor containment 
structures. Several models have been proposed, and they provide insight into the physics of the aircraft 
impact problem (e.g., Riera 1980; Gardner 1984; Zimmerman 1981). For example, in one model, the 
aircraft was viewed as a low-density structure with a crush strength of <j c and an average density of p 
(mass per unit axial length) which vary along the axis of the aircraft with axial position ( x ). 
Furthermore, it was assumed that as the aircraft impacted the building, the aircraft crushed from front to 
back while the building crushed very little. In this way, the force on the building was composed of a 
contribution from the strength of the aircraft and a contribution from the change in momentum of the 
material entering the crush zone. With these assumptions, a simple expression for the force exerted on a 
rigid structure during an impact by the aircraft was 



F(t) = <r c (x) + p(x) 



dx 

~dt 



(1) 



NISTNCSTAR 1-2B, WTC Investigation 369 



Chapter 10 

where x = x(t) is the position of the crush zone in the aircraft as a function of time. The first term in 
Equation (1) is the contribution from the strength of the aircraft, and the second is the contribution from 
the momentum change of the aircraft and its contents. In the Riera approach, the second term in Equation 
(1) was multiplied by a constant factor to obtain better correlation with test data. These types of models 
have often been used to show that the contribution of aircraft strength is on average much smaller than the 
contribution from momentum change. 

Considering the aircraft impact conditions for the WTC towers, the Boeing 767-200ER has a length of 
approximately 259 ft. The empty weight of the Boeing 767 is about 180,000 lb, with a weight at impact 
of approximately 280,000 lb. Much of the added weight of 100,000 lb at take-off was fuel (about 60,000 
lb), with passengers and cargo making up most of the rest. At 500 mph, an impact on a rigid structure 
takes place over a period of about 350 ms (aircraft length divided by aircraft speed). The global impact 
analyses described in Chapter 9 show that this is a good approximation to the duration over which the 
majority of the aircraft momentum was transferred to the towers (although much of the impact and 
momentum transfer is at depth inside the tower). The average force associated with the impact is the 
change in momentum divided by the pulse length or: 

M*V , 

F= = 1.8xl0 7 /bs (2) 

At v ' 

In contrast, an upper bound on the average strength of the aircraft structure can be obtained by estimating 
the yield strength of the structural mass. The empty fuselage weight was assumed to be 100,000 lb, with 
the remaining 80,000 lb of the empty weight consisting of engines (18,000 lb), wings, and landing gear. 
In addition, it was assumed that 20 percent or 20,000 lb of the fuselage weight was structural aluminum 
with non-structural components such as electronics, air conditioning, power units, ductwork, electronic 
wiring, actuator motors, insulation, hydraulics, galley, and lavatories, etc. It was also assumed that this 
structural mass was uniformly distributed along the length of the aircraft. The uniform structural weight 
per unit length was, thus, 77 lb/ft. The aluminum has a density of 175 lb/ft 3 , resulting in an average cross- 
sectional area of 63 in 2 . Assuming the entire aluminum cross section was at the material yield strength of 
50 ksi resulted in an average crush strength of 3.1x10 s lb, or about 17 percent of the inertial force, from 
momentum change. Thus, from this simplified analysis, the contribution of the strength and fracture 
effects should be considerably smaller than inertial effects on the overall impact load. 

This simplified analysis confirms the conclusions of previous studies that the momentum transfer is more 
important than the strength of the impacting aircraft in determining the overall loads on the building. 
However, the finite element simulations show that the strength and geometry of some components of the 
aircraft play an important role in the subsequent damage to the building. The engines are an important 
example, since they are heavy and strong and can penetrate deep into the building. The finite element 
simulations show that large fragments from heavy structural members such as the wing main spar, 
wheels, and similar heavy components penetrate the building and change the nature of the damage. 

10.3 ROTATIONAL AND TRANSLATIONAL ENGINE KINETIC ENERGY 

The model of the aircraft engines used in the global finite element impact analyses did not account for the 
rotational velocity of the spinning components. The additional model complexity required to support the 
inertial loads for the spinning components (as well as clearances, bearing surfaces, etc.) was not believed 



370 NISTNCSTAR 1-2B, WTC Investigation 



Supporting Simplified Analyses 



to be warranted for evaluation of the impact damage. As a result, the rotational kinetic energy of the 
spinning aircraft engine components was not included in the impact analyses. 

An analysis was performed to quantify the potential effect of the rotational kinetic energy of the engine 
turbine and compressor components on the aircraft impact response. This was investigated using simple 
energy and momentum principals. The inertial properties of the engine components were calculated from 
the finite element engine model using LS-DYNA. The inertial properties were automatically generated 
for each part (material assignment) by LS-DYNA as part of the model initialization. 

The internal engine components for which some portion would be rotating are shown in Figure 10-1. 
The engine components in the model were separated to estimate the rotating engine mass and obtain an 
estimated percentage of the engine that is rotating. The approximation used was that 3,500 lb, or 
37 percent of the total engine mass, was rotating. 





Figure 10-1. Internal components of the aircraft engine model. 

Both the kinetic energy from translation and rotation were calculated based on an assumed impact speed 
of 500 mph and an engine rotational speed of 5,000 RPM. The results are shown in Table 10-1. The 
initial rotational kinetic energy of the engine components was approximately 2.2><10 7 lb-ft, which was 
approximately one-quarter of the translational kinetic energy for the engine or less than one percent of the 
initial kinetic energy for the aircraft impact. 





Table 10-1. 


Kinetic energy partitioning in aircraft. 




Component 


Weight 


Moment of Inertia 
about Engine Axis 


Translational 
Kinetic Energy 


Rotational Kinetic 
Energy 


Total Engine 


9400 lb 


3.3x104 lb-ft 2 


8.1x107 lb-ft 


3.5x106 lb-ft 


Rotor Components 


3500 lb 


5.2x103 lb-ft 2 


3.0x107 lb-ft 


2.2x107 lb-ft 


Entire Aircraft 


290,000 lb 




2.4x109 lb-ft 





Since the rotational kinetic energy is a nontrivial percentage of the total engine kinetic energy, additional 
consideration was given to the rotational effects on the engine impact response. The concern is the 
coupling of the rotational energy into the tower and the potential of the rotating engine parts to stabilize 
the engine trajectory during the impact. However, as the engine impacted and penetrated the exterior wall 
of the tower, a significant amount of debris from the tower steel, glass, aluminum trim, and concrete floor 
slab was forced into the engine. It is expected that this would produce significant internal damage and 
lock up the rotating engine components in a relatively short duration. The forces produced by the internal 
angular deceleration of the rotating engine components would either be reacted by the engine support 
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structures or cause the entire engine to roll about its major rotational axis. Using conservation of angular 
momentum for the engine, it was estimated that the rotational energy after the engine locked up was about 
3.5><10 6 lb-ft, or 16 percent of the initial rotational kinetic energy. Therefore, the bulk of the initial engine 
rotational energy would most likely be dissipated by deformation of the internal engine components rather 
than damage to the tower. As a result, the potential for the rotational kinetic energy to significantly 
change the magnitude of the tower impact damage is small. 



10.4 



EFFECTS OF STATIC STRUCTURAL PRELOAD 



The global impact analyses required the use of very large and complex tower models of approximately 
1.5 million nodes. The global WTC tower impact analysis run times were on the order of two weeks on 
twelve processors for a 0.65 second impact event. Initialization of a static preload in these tower models 
was difficult and computationally intensive. As a result, an analysis was performed to determine the 
importance of the preload effect on the tower impact response and damage. The internal energy 
associated with the elastic service loads was small compared to the material internal energy capacity that 
was utilized under impact response and damage. Therefore, the static preload was not expected to have a 
significant influence on the dynamic impact response and deformation. 

A wide flange core column was used for this analysis, with the applied loading and boundary conditions 
as shown in Figure 10-2. The column dimensions were a 16-in. flange width and 16-in. length between 
outer flange faces, with a flange thickness of 1.72 in. and a web thickness of 1.07 in. These dimensions 
were typical for the range of wide flange core columns in the impact zones of the WTC towers. A static 
compressive axial load, P sta tio was applied to the upper end of the column, as shown in Figure 10-2(a). 
An impulsive load per unit area, Pdynamic, was centered between floors at the middle of the column over a 
distance of 60 in. The pressure was applied perpendicular to the flange. The lateral constraints allowed 
the ends of the horizontal beam sections connected to the column to translate only in the direction of the 
static load, as shown in Figure 10-2(b). 



Horizontal booms 

constrained in x. y 

directions 



r-u 



Floor 
Height 



Roller Boundary 
^-translation only 
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(a) Loading detail 



u 

(b) Lateral constraint detail 



Figure 10-2. Problem geometry for analysis of axial preload effects. 
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Three cases were analyzed to study the effects of the axial preload: 

Case 1: A static axial load was applied with no impulsive pressure (Pdynamic = 0.0). The column 
collapsed when the axial force reached 3,360,000 lbf. 

Case 2: An impulsive load was applied with Pdynamic = 2,538 psi and without a static preload. The 

pressure pulse had a characteristic time of 4 ms, resulting in a total impulse of I = 10.2 psi-s. 
A static axial load was subsequently applied and increased until collapse. The column 
collapsed when axial load reached 2,300,000 lbf. 

Case 3: An axial preload was applied until stresses in the column equaled 40 percent of the column's 
yield strength. Then, the lateral impulsive load used in Case 2 was applied. Finally, the axial 
load was increased until collapse. Collapse occurred when axial load reached 2,240,000 lbf. 

Comparing Case 1 with Case 2, it was found that the damage produced by the dynamic load reduced the 
amount of axial force the column was able to support by about 32 percent. In Case 3, with an axial 
preload, the column required 2.6 percent less axial force for collapse than in Case 2 with no axial preload. 
This result suggests that the time at which the axial load was applied relative to the application of the 
dynamic load had little effect on the amount of axial force required to collapse the column. Figure 10-3 
shows a comparison of the damage and plastic strains for Case 2 and Case 3, immediately after the 
application of the dynamic load. The plastic strains for these two cases were very similar. Therefore, the 
static preload has little effect on the damage developed by a lateral impulsive load. 

These comparisons show that the initialization of a static preload in the columns of the towers would not 
have a significant effect on the tower impact response, damage, or residual strength. 
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(a) Case 4, no static load (b) Case 3, static load 

Figure 10-3. Dynamic column response with and without axial preload. 
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10.5 SIMPLIFIED ANALYSIS OF CORE COLUMN DAMAGE 

In order to better understand the importance of the initial impact and breakup of the aircraft on damage to 
the WTC tower core structures, a model of resultant impact damage from an aircraft debris cloud on core 
columns was developed. This model was based on simple momentum and energy conservation concepts, 
similar to the simplified global analysis discussed in Section 10.2. The fundamental premise of the 
analysis is that the aircraft debris has negligible strength on impact, and the load on the building columns 
is due to momentum transfer. As a consequence, the load on the columns is determined by the debris 
density and impact kinematics. The impact pressure delivered to the column is determined by momentum 
change. The pulse duration is determined by the debris cloud speed and size. This analysis was used to 
study the effects of the dispersed fuel cloud on core column response. 

The critical loads required to fail a core column were determined based on a combination of the peak 
impact pressure and impulse intensity. A good description of this characterization of critical loads as 
Pressure-Impulse or "P-I" damage curves is given by Abrahamson and Lindberg [1972]. The concept that 
broad definitions of critical loads can be made based on peak pressure or load and impulse intensity or 
total impulse is well established. The P-I damage analysis method has been commonly used for 
evaluation of blast damage to structures and is appropriate for structures where the impacting body has a 
low strength. This P-I damage curve approach provides a useful way to make a quick evaluation of the 
damage potential of aircraft debris cloud on core columns. 

10.5.1 Development of P-I Curves 

An example of the process used to develop the P-I curves for a structure is demonstrated in Figure 10-4. 
A series of simple finite element calculations with various pressure and impulse combinations are 
performed to bracket the conditions required to initiate damage or failure of the structure. Each 
calculation is represented by a circle in the figure. The load conditions that resulted in no damage or 
incipient damage are represented by a solid blue circle. Load combinations resulting in failure of the 
structure are represented with a solid red circle. Intermediate levels of damage are displayed with open 
black circles. The results from the simulations are then fit with two P-I curves, one representing the onset 
of significant damage and the other representing complete failure. The P-I curves have a form that is 
similar in shape to a hyperbolic curve. The load combinations which fall above and to the right of the red 
P-I curve in Figure 10-4 produce structural failure. Similarly, load combinations which fall below and to 
the left of the blue P-I curve in Figure 10^4 produce little or no damage. 

Individual P-I curves correspond to a specific magnitude of damage. In this example, separate P-I curves 
are generated for the onset of damage and for failure. A comparison of the two P-I curves in Figure 10-4 
shows that at the limit of a long duration loading, the P-I curves are representing damage and failure are 
equal for this structure (lower right portion of the P-I curves with a low pressure magnitude and large 
impulse intensities). This results from unstable behavior of this structure under quasi-static pressure 
loads. As a result, the quasi-static pressure magnitude that is sufficient to produce yielding will also 
result in a complete rupture. The knee of the P-I curves in this example corresponds to an effective load 
duration between one and ten milliseconds, which corresponds to the characteristic response times of the 
structure (natural frequencies of the important damage modes such as column buckling modes, etc.). 
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Figure 10-4. Development of damage and failure P-l curves for a structure. 

10.5.2 Core Column Damage Analyses 

Using the procedure outlined in the previous section, P-I curves were developed for core column failure. 
The column dimensions and boundary conditions were identical to those used in Section 10.4 and shown 
in Figure 10-2, but without an axial preload. Three different P-I curves were developed for three heights 
of dynamic pressure loading, as shown in Figure 10-5. Pressure and impulse combinations were applied 
over 30 in., 60 in., and 120 in. spans on the column. The P-I curves were then used to determine the 
characteristics of a fuel cloud impact that would be sufficient to fail a core column. 

The core column response was analyzed for the aircraft fuel impact loads at four positions along the wing. 
The impact loading of the fuel was considered separately in this example since it represents a significant 
fraction of the total mass in the wings. Initially, the effect of a direct impact by a section of the fuel on 
the core column (no dispersion of the fuel cloud) was considered. To calculate the impact pressure and 
impulse, the fuel tank dimensions at each location were approximated by a rectangle. Near the root of the 
wing the approximate tank dimensions were 43 in. tall by 161 in. across. The position nearest the wing 
tip had approximate tank dimensions of 14 in. tall by 50 in. across. The other two positions had relatively 
uniform spacing and had dimensions of 18 in. tall by 91 in. across and 22 in. tall by 132 in. across, 
respectively. The peak pressure for the fuel impact was found from the momentum change, which is 
governed by the fuel density and impact velocity. The duration of the load was calculated using the 
dimension across the fuel tank divided by the impact velocity. The pressure multiplied by the duration 
gives the impulse intensity for this impact scenario. 
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The impact loads for the four different positions on the fuel tank were calculated for three impact speeds, 
537, 430, and 268 mph. The corresponding positions of the impact loads in pressure and impulse are 
shown in Figure 10-5. The comparison of the impact loads to the P-I curve for the 30 in. load span 
(green curve) is most appropriate since the loaded area in this example varies from 43 in. at the wing root 
to 14 in. at the wing tip. These data show that a solid section of fuel directly impacting the core columns 
at 430 mph or above would cause failure. However, the loads produced by a lower impact speed of 
268 mph would not be sufficient to fail this core column. 
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Figure 10-5. Core column P-l failure curves and concentrated fuel impact loading. 

The fuel impact loads used in the above example were an upper bound on the damage that could be 
produced by the fuel at the core since the fuel was still assumed to be in a solid rectangular section at full 
density. In reality, the impact with the exterior wall broke up the wing structures, and the fuel cloud 
spreads out to some extent prior to reaching the core columns. The impact with the internal building 
contents also reduced the severity of the fuel impact at the core. 

A first approximation of the effects of the fuel cloud dispersion on the core column loading was 
performed by assuming that the fuel cloud spread vertically to fill the full floor height prior to reaching 
the core. A uniform fuel cloud density was used that was scaled by the ratio of the initial fuel section 
height to the full floor height. No spreading in the longitudinal or lateral direction was assumed. The 
spreading of the debris cloud, thus, reduced the peak pressure without changing the load duration. These 
positions of the impact loads in pressure and impulse are shown in Figure 10-6. The comparison of the 
impact loads to the P-I curve for the 120 in. load span (black curve) was most appropriate since the loaded 
area in this example was uniform over the full length of the column. The comparison shows that an 
expanded fuel cloud did not produce sufficient loading to fail a core column with the exception of points 
close to the wing root at the highest impact speeds considered. A more realistic fuel cloud dispersion 
would have included lateral and longitudinal spreading, as well as removal of some fuel from the cloud as 
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a result of impact with the outer wall and building contents. Under these conditions, it is expected that 
the fuel cloud alone would not be sufficient to fail core columns. 
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Figure 10-6. P-l curves and expanded fuel cloud impact loading. 
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11.1 INTRODUCTION 

The objective of this project was to analyze the aircraft impacts into the World Trade Center (WTC) 
towers to provide the following: (1) estimates of probable damage to structural systems, including the 
exterior walls, floor systems, and interior core columns; (2) estimates of the aircraft fuel dispersal during 
the impact; and (3) estimates of debris damage to the interior tower contents, including partitions and 
workstations to be used for estimating damage to fire proofing and to the mechanical and architectural 
systems inside the towers. This project, thus, established the initial conditions for the fire dynamics 
modeling in Project 5 and the thermal-structural response and collapse initiation analysis in Project 6 of 
the National Institute of Standards and Technology (NIST) Investigation. 

The WTC aircraft impact analysis is very complex with large scale fracture and fragmentation of both 
tower and aircraft structures, nonlinear rate-dependent material behaviors, and the fluid- structure 
interaction of the aircraft fuel. The analyses of the aircraft impacts performed for this investigation are 
believed to be the highest-fidelity simulations ever performed for this impact behavior using state-of-the 
art analysis methodologies. Wherever possible, the models were validated against observables or 
supporting test data developed by the WTC Investigation. 

The impact analyses were conducted at various levels of complexity including: (1) the component level, 
(2) the subassembly level, and (3) the global level to estimate the probable damage to the towers due to 
aircraft impact. Analysis of uncertainties using the component, subassembly, and global analyses were 
also conducted to assess the effect of uncertainties associated with various parameters on the damage 
estimates and to identify the most influential parameters that affect the impact response. 

11.2 MATERIAL CONSTITUTIVE MODELING 

An important requirement for high fidelity simulation of the aircraft impact damage was the development 
of constitutive models that represented the actual behavior of the WTC towers and aircraft under the 
dynamic impact conditions. The materials that were considered included: (1) the several grades of steel 
used in the columns, spandrels, and floor trusses, and beams of the WTC towers, (2) the concrete floor 
slabs, (3) the various aluminum and titanium alloys used in the aircraft, and (4) the nonstructural contents 
of the towers. These materials exhibit significant nonlinear, rate-dependent deformation and failure 
behavior that need to be represented in the constitutive relationship. 

The primary constitutive model that was used for the tower steels was the Piecewise Linear Plasticity 
model in LS-DYNA. This model is sufficient to model the nonlinear rate-dependent deformation and 
failure of the steel structures. The constitutive model parameters for each grade of steel were based on 
engineering stress-strain data obtained by Project 3 of the NIST Investigation. Finite element analyses of 
the test specimens were conducted with a fine and a medium mesh to capture the nonlinear material 
behavior up to failure. The finite element analysis provided a validation that the constitutive model 
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parameters were defined accurately and that the model could reproduce the measured response for the test 
conditions. Strain-rate effects on the steel yield strength were included in the constitutive model for 
tower steels with the Cowper and Symonds rate effect model. 

The LS-DYNA material Type 16 (pseudo-tensor concrete model) was selected for modeling the concrete 
floor slabs due to its ability to accurately model the low confinement damage and softening behavior of 
concrete. The model uses two pressure-dependant yield functions and a damage-dependent function to 
migrate between curves. This allows for implementation of tensile failure and damage scaling, which are 
more dominant material behaviors at low confinement. The pseudo-tensor model also accounts for the 
high strain-rate sensitivity of concrete. 

The constitutive and failure properties for the aircraft materials were developed from data available in the 
open literature. Complete engineering stress-strain curves were obtained for various 2024 and 7075 
aluminum alloys that are commonly used in the construction of the Boeing 767 airframe structures. 
These curves were digitized for the various aluminum alloys. Representative stress-strain curves were 
then converted into true stress and true strain and used to develop tabular curves for the constitutive 
models. No rate sensitivity of the aircraft materials was considered. 

The primary influence of the nonstructural components on the impact behavior was their inertial 
contribution. The effects of their strength were small. As a result, relatively simple approximations of 
their constitutive behavior were used. Typically a bilinear elastic-plastic constitutive model was applied 
for these materials to allow for efficient modeling of deformation and subsequent erosion from the 
calculations as their distortions became large. 

11.3 TOWER MODEL DEVELOPMENT 

One of the significant challenges in developing the tower and aircraft models for the global impact 
analyses was to minimize the model size while keeping sufficient fidelity in the impact zone to capture 
the deformations and damage distributions. The limitation was that the combined aircraft and tower 
models should not exceed approximately 2.3 million nodes. These were distributed between the global 
WTC tower model and the aircraft so that the tower model would be about 1.5 million nodes and the 
aircraft about 0.8 million nodes. The approach used to meet this objective was to develop models for the 
various tower components at different levels of refinement. Components in the path of the impact and 
debris field were meshed with a higher resolution to capture the local impact damage and failure, while 
components outside the impact zone were meshed more coarsely to primarily capture their structural 
stiffness and inertial properties. 

The towers were modeled primarily with shell elements, with the exception of the exterior wall bolted 
connections (beam and brick elements) and the floor truss diagonals (beam elements). The WTC 1 model 
extended between floors 92 and 100, while the WTC 2 model extended between floors 85 and 77. The 
global impact models of the WTC towers included the following components: 

• Core columns and floors: Core columns were modeled using shell elements with two mesh 
densities, a refined density in the direct impact area and a coarser far field density elsewhere. 
The spliced column connections were included in the model with proper failure criteria. The 
floors within the core were modeled using shell elements representing the floor slabs and 
beams. 
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Exterior walls: The exterior columns and spandrels were modeled using shell elements with 
two mesh densities, a refined density in the immediate impact zone and a coarser far field 
density elsewhere. For the bolted connections between exterior panels in the refined mesh 
areas, brick elements were used to model the butt plates, and beams elements were used for 
the bolts. 

• Truss floor: In the direct impact area, the floor model included shell elements for the 
combined floor slab and metal decking and for the upper and lower chords of the trusses. 
Beam elements were used for the truss diagonals. In the far field floor segments, simplified 
shell element representations were used for the floor slab and trusses. The dampers were not 
included in the models due to their low mass and strength, particularly in the transverse 
direction. 

• Interior building contents: The interior nonstructural contents of the towers were modeled 
explicitly. These included the partitions and workstations, which were modeled with shell 
elements in the path of the aircraft debris. The live load mass was distributed between the 
partitions and cubicle workstations. 

11.4 AIRCRAFT MODEL DEVELOPMENT 

The finite element model for the Boeing 767-200ER aircraft was constructed through a three-step process: 
(1) data collection, (2) data interpretation and engineering analysis, and (3) meshing of the structure. The 
focus of this effort was on gathering sufficient structural data and including adequate detail in the aircraft 
model so that the mass and strength distribution of the aircraft and contents were properly captured for 
implementation in the impact analyses. Structural data were collected for the Boeing 767-200ER aircraft 
from (1) documentary aircraft structural information and (2) data from measurements on Boeing 767 
aircraft. 

The airframe model contained most of the significant structural components in the aircraft. The models 
of the fuselage, empennage, and wing structures were developed completely using shell elements. 
Models for the landing gear and engines were primarily developed using shell elements, but contained 
some brick elements as well. The typical element dimensions were between one and two in. for small 
components, such as spar or rib flanges, and three to four in. for large parts, such as the wing or fuselage 
skin. 

Special emphasis was placed on modeling the aircraft engines due to their potential to produce significant 
damage to the tower components. The engine model was developed primarily with shell elements. The 
objective was to develop a mesh with typical element dimensions between one and two inches. However, 
smaller element dimensions were required at many locations to capture details of the engine geometry. 
Brick elements were used for some of the thicker hubs and the roots of the compressor blades. 

A detailed analysis was carried out to estimate the fuel distribution in the aircraft wings at the time of 
impact. 
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11.5 COMPONENT LEVEL ANALYSES 

The primary objectives of the component modeling were to (1) develop understanding of the interactive 
failure phenomenon of the aircraft and tower components and (2) develop the simulation techniques 
required for the global analysis of the aircraft impacts into the WTC towers. The approach taken for 
component modeling was to start with finely meshed, brick and shell element models of key components 
of the tower and aircraft structures and progress to relatively coarsely meshed beam and shell element 
representations used in the global models. In addition to determining the optimal element size and type 
for global modeling, other key technical areas were addressed in the component modeling, including 
material constitutive modeling, treatment of connections, and modeling of aircraft fuel. Examples of the 
component impact analyses included: 

Impact of a segment of an aircraft wing with an exterior column. 

Impact of a segment of an aircraft wing with a core column. 

Detailed and simplified modeling of exterior panel bolted connection under impact loading 
and modeling of the bolted spandrel connection. 

Impact of an aircraft engine with exterior wall panels. 

Impact of a simplified plow type impactor with truss floor assembly. 

Impact of an empty wing segment with exterior wall panels. 

Impact of fuel-filled wing segment with exterior wall panels. 

The following results were obtained from the component impact analyses: 

• A 500 mph engine impact against an exterior wall panel resulted in a penetration of the 
exterior wall and failure of impacted exterior columns. If the engine did not impact a floor 
slab, the majority of the engine core would remain intact through the exterior wall penetration 
with a reduction in speed between 10 percent and 20 percent. The residual velocity and mass 
of the engine after penetration of the exterior wall was sufficient to fail a core column in a 
direct impact condition. Interaction with additional interior building contents prior to impact 
or a misaligned impact against the core column could change this result. 

• A normal impact of the exterior wall by an empty wing segment from approximately mid- 
span of the wing produced significant damage to the exterior columns but not complete 
failure. Impact of the same wing section, but filled with fuel, resulted in extensive damage to 
the external panels of the tower, including complete failure of the exterior columns. The 
resulting debris propagating into the building maintained the majority of its initial momentum 
prior to impact. 

• Three different numerical techniques were investigated for modeling impact effects and 
dispersion of fuel: (1) standard Lagrangian finite element analysis with erosion, (2) Smoothed 
Particle Hydrodynamics (SPH) analysis, and (3) Arbitrary-Lagrangian-Eulerian (ALE) 
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analysis. Of these approaches, use of the SPH offered the best viable option due to its 
computational efficiency. 

11.6 SUBASSEMBLY IMPACT ANALYSES 

The subassembly analyses were considered as a transition between the component level analyses and the 
global impact analyses. With the subassembly analyses, more complex structural behavior not captured 
in the component analyses could be investigated with significantly shorter run times than required for the 
global analyses. The subassembly analyses were primarily used to investigate different modeling 
techniques and associated model size, run times, numerical stability, and impact response. The final 
subassembly model used structural components from the impact zone on the north face of WTC 1, and 
included exterior panels, core framing, truss floor structures, and interior contents (workstations). The 
subassembly model was impacted by an aircraft engine and by a segment of a fuel-filled wing. 

The subassembly model was used to investigate the effect of a number of modeling parameters on the 
response and damage estimates. For the engine impact simulations, these parameters included the 
strength of the building nonstructural contents and the concrete slab strength. For the wing impact 
simulations, the effect of the ductility of the exterior column weldment on the impact response was 
investigated. The results of these parametric studies indicated the following: 

• The deceleration profile of the impacting engine indicated that the response of the 
nonstructural building contents was dominated by the mass of the workstations, rather than 
by their strength. 

• Varying the strength of the floor concrete slab from 4 ksi to 3 ksi did not result in significant 
change in the impact response. It appears that the mass of the concrete slab had a greater 
effect on the engine deceleration and damage to the floor than did the concrete strength. 

• Varying the ductility of the weld zone in the exterior columns from 8 percent to 1 percent did 
not result in any noticeable difference in the damage pattern or the energy absorbed by the 
exterior panels, indicating that the weld ductility had a negligible effect on the impact 
response. 

11.7 ANALYSIS OF AIRCRAFT IMPACT CONDITIONS 

The initial aircraft impact conditions included aircraft speed, horizontal and vertical angles of incidence, 
roll angle of each aircraft, and the location of nose impact with each tower. Three different methods were 
applied to determine the impact conditions for the two aircraft that impacted the towers. The first method 
used a comparison of videos from different positions to calculate the three-dimensional trajectory of the 
aircraft. The second method used the relative frame-by-frame motion in a single video scaled to the 
length of the aircraft in the video to calculate the impact speed. Finally, analysis of the impact damage on 
the face of each tower was used to refine the relative impact orientation and trajectory. This was done by 
matching the projected impact points of the wings, fuselage, engines, and vertical stabilizer onto the 
exterior wall of each tower to the observed damage pattern. 
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The analysis indicated that the impact speed on WTC 2 (542 ± 24mph) was about 23 percent larger than 
the impact speed on WTC 1 (443 ± 30 mph). The aircraft that impacted the north tower was steeper than 
the aircraft that impacted the south tower (vertical approach angle of 10.6° ± 3° for WTC 1 versus 6° ± 2° 
for WTC 2). The roll angle for WTC 1 was 25° ± 2° while that for WTC 2 was 38° ± 2°. 

11.8 UNCERTAINTY ANALYSES 

The objectives of the uncertainty analyses were to assess the effect of uncertainties associated with the 
aircraft and WTC towers parameters on the level of damage to the towers after impact and to determine 
the most influential modeling parameters that affect the damage estimates. Uncertainty arises in these 
analyses from the following key parameters: (1) aircraft impact parameters (speed, horizontal and vertical 
angles of incidence, orientation, and location of impact), (2) material properties and failure criteria for the 
towers and the aircraft, (3) aircraft parameters (mass and stiffness properties, and jet fuel distribution), 

(4) tower parameters (structural strength and mass distribution, connections behavior), and 

(5) nonstructural building contents that may share in absorbing energy imparted by the aircraft impact. 

Another important source of uncertainty is the inaccuracy associated with mathematical or numerical 
models. These uncertainties, also known as modeling errors, are deterministic in nature, but are often 
treated as random variables to characterize the effects of the analysis methodologies on the calculated 
response. All of these variables did not necessarily have a significant effect on the estimated impact 
damage to the WTC towers. 

Because of the complexity of the problem and the limited number of parameters that could be varied in 
the global analyses, it was necessary to down-select a refined list of uncertainty parameters from all of the 
possible parameters. Therefore, parameter screening was conducted using design of experiments 
methodology. Screening was first conducted at component and subassembly levels using orthogonal 
factorial design techniques in order to identify the most influential parameters and reduce the number of 
parameters to a more manageable number for the global impact analyses. The sensitivity analyses 
included engine impacts against core columns, wing section impacts against exterior panels, and engine- 
impact subassembly analyses. 

Based on the three sensitivity analyses, the set of influential modeling parameters was reduced. The 
following parameters were selected for variation in the global impact analyses: 

Impact speed. 

Vertical approach angle of the aircraft. 

Lateral approach angle of the aircraft. 

Total aircraft weight. 

Aircraft materials failure strain. 

Tower materials failure strain. 

Building contents weight. 
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The impact speed and vertical approach angle were selected as significant parameters in the global 
analyses due to their importance in the component and subassembly impact response. The vertical 
approach angle played a primary role in the magnitude of the vertical impact loads on the truss floor 
structures. Lateral approach angle was also selected, as this parameter dictated to a large extent where 
aircraft debris traveled and what part of the core was affected by this debris. 

Uncertainties in the strength and ductility (failure strain) of materials had a coupled effect. An increase in 
yield strength or an increase in failure strain resulted in an increase in the energy absorbing capacity of 
the structure. The yield strength of materials was typically known more accurately than the failure 
strains. This was particularly true within the finite element analyses, where the value of the failure strain 
needed to be assigned based on the model resolution and failure criteria used. As a result, only the 
uncertainties in the material failure strain were used as a material uncertainty parameter for both the 
aircraft and the towers in the subsequent global impact analyses. 

The uncertainty in the weights associated with building contents (corresponding to service live loads) was 
found to be of secondary importance in the engine-subassembly impact analysis. In the global impact 
analyses, the live load contents were expected to play a more significant role in confining the fuel and 
debris dispersion. In addition, the partition walls were significant for controlling the subsequent spread of 
fire through the towers. As a result, the uncertainty in building contents weight was included as an 
uncertainty parameter in the global analyses. 

11.9 GLOBAL IMPACT ANALYSES 

The primary objective of this project was to determine the condition of the WTC towers immediately 
following the aircraft impacts. This assessment included the estimation of the structural damage that 
degraded their strength and the condition and position of nonstructural contents such as partitions, 
workstations, aircraft fuel, and other debris that influenced the behavior of the subsequent fires in the 
towers. The global impact analyses were the primary method by which the damage to the towers was 
estimated. The global analyses included, for each tower, a "base case" based on a best estimate of all 
input parameters. They also provided more and less severe damage estimates based on variations of the 
most influential parameters. These more and less severe damage scenarios provided a range of damage 
estimates for the towers due to aircraft impact. 

11.9.1 WTC 1 Base Case Global Impact Analysis 

The base case impact analysis was performed for a 0.715 s duration following initial impact of the aircraft 
nose with the north exterior wall. The analysis was performed on a computer cluster using twelve 
2.8 GHz Intel Xeon processors, each on a separate node of the cluster. The run time for this analysis was 
approximately two weeks. The residual kinetic energy of the airframe components at the termination of a 
global impact simulation was typically less than one percent of the initial kinetic energy at impact. 

The aircraft impact response was dominated by the impact, penetration, and fragmentation of the airframe 
structures. The entire aircraft fully penetrated the tower at approximately 0.25 s. The fuselage structures 
were severely damaged both from the penetration through the exterior columns and the penetration of the 
96th floor slab that sliced the fuselage structures in half. The downward trajectory of the aircraft 
structures caused the airframe to collapse against the floor, and the subsequent debris motion was 
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redirected inward along a more horizontal trajectory parallel to the floor. The downward trajectory of the 
aircraft structures transferred sufficient vertical load that the truss floor structures on the 95th and 
96th floors collapsed in the impact zone. 

The wing structures were completely fragmented by the exterior wall. The aircraft fuel cloud began to 
spread out after impact but remained relatively dense until the leading edge of the fuel reached the tower 
core. The aircraft fuel and debris cloud eventually penetrated most of the distance through the core before 
their motion was halted. 

Exterior Wall Damage 

The calculated and observed damage in the impact damage zone were in good agreement. The agreement 
in the position and shape of the impact damage provided partial validation of the modeled geometry of the 
aircraft, including the aircraft orientation, trajectory, and flight distortions of the wings. 

The calculated magnitude and mode of impact damage on the exterior wall also agreed well with the 
observed damage. The exterior wall completely failed in the regions of the fuselage, engine, and fuel- 
filled wing section impacts. Damage to the exterior wall was observed out to the wing tips, but the 
exterior columns were not completely failed in the outer wing and vertical stabilizer impact regions. 
Failure of the exterior columns occurred both at the bolted connections between column ends and at 
various locations in the column, depending on the local severity of the impact load and the proximity of 
the bolted connection to the impact. The agreement of both the mode and magnitude of the impact 
damage served to partially validate the constitutive and damage modeling of the aircraft and exterior wall 
of the tower. 

Core Damage 

The core had significant damage in the region close to the impact point. The columns in line with the 
aircraft fuselage failed on the impact side, and several of the core beams were also severely damaged or 
failed in the impact zone. The calculated damage to the core of WTC 1 consisted of three severed 
columns and four heavily damaged columns. 

Truss Floor Damage 

The trusses experienced significant damage in the impact zone. The calculated impact response produced 
severe damage to the truss structures in the primary impact path of the fuselage. The truss structures were 
severely damaged from the exterior wall to the core. The truss floor system on floors 94 through 96 were 
damaged and sagged downward as a result of the impact loading. 

The truss floor and associated slab was heavily damaged or completely destroyed in the impact zone 
ahead of the fuselage. The zone of heavy damage to the floor slab extended beyond the truss floor and 
was approximately one-third of the distance through the core on floors 94 through 96. 
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Fuel and Debris Distribution 

The bulk of the fuel and aircraft debris was deposited on floors 93 through 97, with the greatest 
concentration on floor 94. The calculated debris cloud included 17,400 lb of debris and 6,700 lb of 
aircraft fuel outside of the tower at the end of the impact analysis, either rebounding from the impact face 
(north wall) or passing through the tower (south wall). This amount might have been over-estimated in 
the calculation since the exterior walls were not modeled with windows that would contain the fuel cloud 
and other small debris inside the towers. In addition, the impact behavior of the aircraft fuel cloud did not 
include the ability to stick to, or wet, interior components. Rather, the aircraft fuel SPH particles tended 
to bounce off of internal structures. 

11.9.2 WTC 2 Base Case Global Impact Analysis 

The WTC 2 base case impact analysis was performed for a 0.62 s duration following initial impact of the 
aircraft nose with the south exterior wall. Full penetration of the aircraft into the tower was completed at 
0.2 s after impact. The aircraft impact response was very similar to that of the WTC 1 impact and was 
dominated by the penetration and fragmentation of the airframe structures. The fuselage structures were 
severely damaged both from the penetration through the exterior columns and the penetration of the 81st 
floor slab that sliced the fuselage structures in half. The downward trajectory of the aircraft structures 
caused the airframe to collapse against the floor, and the subsequent debris motion was redirected inward 
along a more horizontal trajectory parallel to the floor. The downward trajectory of the aircraft structures 
transferred sufficient vertical load that the truss floor structures on the 80th and 81st floors began to 
collapse in the impact zone by the end of the simulation. 

The aircraft wing structures and fuel tank were fragmented by the impact with the tower exterior. The 
aircraft fuel cloud started to spread out immediately after impact, but the leading edge of the fuel 
remained relatively dense until passing approximately one-third of the lateral distance through the tower 
core (approximately 0.2 s after impact). At 0.3 s after impact, the aircraft fuel cloud had penetrated 
approximately two-thirds the distance through the core and was spreading out. Beyond this time, the 
subsequent motion of the aircraft fragments and fuel debris cloud was noticeably slowed. The spread of 
the fuel and debris cloud was more rapid and extensive in the open truss floor regions than through the 
core as a result of the open volume above the workstations in the truss floor zone. 

Exterior Wall Damage 

The exterior wall completely failed in the regions where the fuselage, engine, and fuel-filled wing section 
impacted the structure. Damage to the exterior wall extended to the wing tips, but the exterior columns 
were not completely failed in the outer wing and vertical stabilizer impact regions. Failure of the exterior 
columns occurred both at the bolted connections between column ends and at various locations in the 
column, depending on the local severity of the impact load and the proximity of the bolted connection to 
the impact. 

The calculated and observed geometry and magnitude of impact damage were in good agreement. That 
served to partially validate the geometry of the aircraft model, including the aircraft orientation, 
trajectory, and flight distortions of the wings. Agreement of both the mode and magnitude of the impact 
damage partially validated the constitutive and damage modeling of the aircraft and exterior wall of the 
tower. 
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Core Damage 

In the base case WTC 2 impact damage analysis, significant damage to the core was calculated in the 
region close to the impact point (the southeast corner of the core). The columns in line with the aircraft 
fuselage failed on the impact side of the tower, and several of the core beams were also severely damaged 
or failed in the impact zone. In some cases, failure of the column splices located on floors 77, 80, and 
83 contributed significantly to the failure of the core columns. This was particularly true for the heavy 
column number 1001 at the southeast corner of the core that failed at three splice locations. A total of 
five columns were severed and four columns were heavily damaged. 

Truss Floor Damage 

The trusses experienced significant damage in the impact zone, with the greatest damage on floor 81. 
This severe damage to the floor 81 truss structure and floor slab was in the primary impact path of the 
fuselage and extended from the exterior wall to the core. The truss floor system on floors 79 and 81 had 
sufficient damage from the impact that truss floor sections sagged downward as a result of the impact. 
Trusses on floor 80 also sustained significant damage as a result of impact. 

The calculated truss floor damage to WTC 2 was less than that of WTC 1, despite the higher aircraft 
impact energy for WTC 2. The larger truss floor damage and deflection in WTC 1 can be explained by 
two factors. The primary factor was that the WTC 1 downward impact trajectory was nearly twice as 
steep as that of the WTC 2 impact. As a result, the steeper impact angle directed more impact energy 
normal to the floor slab. The secondary factor was that the damage to the long span truss floor in the 
WTC 1 impact zone produced larger displacements than the corresponding damage level to the short span 
truss region in WTC 2. 

Fuel and Debris Distribution 

The bulk of the fuel and aircraft debris was deposited in floors 78 through 80, with the largest 
concentration of aircraft debris on floor 80, and the largest concentration of aircraft fuel on floors 79, 81, 
and 82. The calculated debris distribution included 55,800 lb of debris and 10,600 lb of aircraft fuel 
outside of the tower at the end of the impact analysis, either rebounding from the impact face or passing 
through the tower. Similar to WTC 1, the calculated debris outside the tower is believed to be larger than 
is realistic. 

11.9.3 Different Severity Global Impact Analyses 

Additional impact analyses were performed for each tower to provide a range of damage estimates due to 
the uncertainties in the calculated impact response. Two additional global impact analyses were 
performed for each tower in order to evaluate the impact response for a more severe and less severe 
impact scenario. The variations in impact conditions for the different severity global impact analyses 
were primarily obtained from the uncertainty analyses. 

WTC 1 Global Impact Analyses 

The impact speeds were 414 mph and 472 mph in the less severe and more severe impact scenarios, 
respectively. These were the upper and lower bounds obtained from the analysis of aircraft impact 
conditions. The vertical trajectory of the aircraft was also varied from 13.6 degrees in the less severe case 
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to 7.6 degrees for the more severe impact scenario. The reduced vertical trajectory angle resulted in a 
greater impact energy directed inward toward the core. The lateral trajectory was not varied since the 
impact was close to centered on the tower and normal to the north face of WTC 1. A small variation in 
the lateral approach angle would have little effect on the energy of the aircraft debris entering the tower 
and core. 

A comparison of the calculated damage to the north exterior wall of WTC 1 from the three different 
severity impact analyses with the observed damage indicated that the overall agreement with the observed 
damage was good for all three analyses, with the base case global impact analysis providing the best 
match to the observed damage. The calculated damage magnitude was similar in each of the global 
analyses. The small differences in apparent damage were largely due to panels that may have severed 
columns in one case and were removed at the connections in another. In general, the trend was for a 
larger opening produced by the less severe impact and a smaller opening in the more severe impact. The 
increase in the opening with reduced severity impact can be explained by the increased tower panel 
material toughness in the less severe scenario, resulting in a transition from severing of columns to the 
failure of connection bolts and panel removal. 

Toward the wing tips, where the columns and spandrels were not completely severed, the more severe 
impact damage analysis calculated larger damage to the exterior wall panels. These columns had the 
largest amount of material with plastic strains above 5 percent. As would be expected, the less severe 
impact damage analysis calculated lower damage to the exterior wall, and the base case analysis 
calculated an intermediate level of damage near the wing tips. 

The core damage was concentrated in the region of the core closest to the fuselage impact point, and there 
was a clear correlation in damage magnitude with the impact severity. A total of one column was 
severed, and two columns were heavily damaged in the less severe impact, compared to three columns 
severed and four columns heavily damaged in the base case impact analysis and six columns severed and 
three columns heavily damaged in the more severe WTC 1 impact analysis. A strong correlation in core 
damage with impact severity was expected. All of the parameter variations in the three increasing 
severity analyses were selected to produce an increase in core damage. A higher impact speed and a 
shallower impact angle were selected to direct more energy toward the core. The aircraft mass was 
increased and given a higher material toughness. The mass of the tower contents was reduced, and the 
tower materials were given a lower toughness. All of these variations contributed toward the increased 
core damage with impact severity. 

All of the WTC 1 global impact analyses resulted in a similar amount of damage to the truss floor and 
floor slab. Competing parameters canceled their respective effects on the floor damage. Variations in 
many of the parameters should have caused greater damage to the floor system in the more severe impact 
analysis. However, the downward impact trajectory angle was also reduced to direct more energy into the 
core. Therefore, less of the impact energy was oriented into the floor system. The combined effects of 
varying these parameters led to very similar damage to the floor system in all of the analyses. 

WTC 2 Global Impact Analyses 

The impact speeds were 521 mph and 570 mph in the less severe and more severe impact scenarios, 
respectively. These speeds were the upper and lower bounds obtained from the analysis of aircraft impact 
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conditions. The vertical trajectory of the aircraft was also varied from 5 (more severe case) to 8 degrees 
(less severe case) in order to direct less energy into the core for the less severe case. 

Calculated damage to the exterior wall and truss floor of WTC 2 showed similar trends to those seen in 
the WTC 1 analyses. Only slight changes in damage to these structures were calculated from the different 
severity impacts. Damage to the exterior was relatively insensitive to the parameter variations. Truss 
floor damage did not vary significantly due to the competing effect of the parameters varied. 

The core damage was concentrated in the core region closest to the impact point. There is also a clear 
correlation in damage magnitude with the impact severity. A total of three columns were severed, and 
two columns were heavily damaged in the less severe impact, compared to five columns severed and four 
columns heavily damaged in the base case analysis and ten columns severed and one columns heavily 
damaged in the more severe WTC 2 impact analysis. 

11.9.4 Comparison with Observables 

The observables available to help validate the global impact analyses included the following: 

• Damage to the building exterior (exterior walls and floors in the immediate vicinity of the 
impact) documented by photographic evidence. 

• Aircraft debris external to the towers (landing gear for WTC 1 and a landing gear and an 
engine for WTC 2) as documented by photographic evidence. 

• Eyewitness accounts from survivors who were inside the towers (blocked or passable 
stairwells). 

Not all of these observables were perfectly matched by the simulations due to the uncertainties in exact 
impact conditions, the imperfect knowledge of the interior tower contents, the chaotic behavior of the 
aircraft break up and subsequent debris motion, and the limitations of the models. In general, however, 
the results of the simulations matched these observables reasonably well. 

11.9.5 Comparison with Previous Studies 

Two studies were previously conducted to estimate the impact damage to the WTC towers. These studies 
were performed by staff at Massachusetts Institute of Technology (MIT) and Weidlinger Associates, Inc. 
(WAI). The MIT study used an energy balance approach to estimate damage to the core columns. The 
WAI study used the FLEX finite element code to calculate the aircraft impact damage to both towers. 

In both towers, the impact damage to the exterior wall calculated in this study closely matched the 
observed damage. The damage profiles in the WAI impact simulations had some noticeable differences. 
The first was that damage predictions from WAI included complete failure of the exterior columns over 
the entire length of the wings and top of the vertical stabilizer. The second was that the WAI predictions 
displayed local rupture of the columns adjacent to the impact point with less influence of the bolted 
connections on panel failure. Neither is consistent with the observed damage. 
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The differences in the damage profiles between the WAI and NIST calculations most likely resulted from 
a variety of differences in the aircraft and tower models. The WAI aircraft model was based on their 
model of a Lockheed C-141B military transport rather than the Boeing 767. The NIST aircraft model also 
contained an order of magnitude more elements and explicit modeling of the fuel. There were also 
differences in the aircraft impact conditions. 

In general, the MIT and WAI studies appeared to predict more damage to the core columns compared to 
the NIST estimates. For WTC 1, MIT estimated that between four to 12 core columns failed. WAI gave 
two estimates for core column failure. The first estimate of 23 core columns failed and five damaged was 
obtained from the FLEX impact analysis. The second estimate of 20 failed columns was the number used 
in their collapse analysis. The NIST base case impact damage of three severed and four heavily damaged 
core columns and less severe estimate of one severed and two heavily damaged core columns fell below 
both the MIT and WAI estimates. The more severe estimate of six severed and three heavily damaged 
columns fell in the middle of the MIT range, but still well below the WAI estimates. 

A similar trend in the predicted damage to the core columns was found in the WTC 2 analysis. MIT 
estimated between seven to 20 columns failed. WAI calculated 14 core columns failed and another 
10 damaged in their FLEX analysis, but reduced the number of failed columns to five for their collapse 
analysis. The NIST base case impact damage of five severed and four heavily damaged columns, as well 
as the more severe estimate of 10 severed and one heavily damaged columns fell in the middle of the 
range predicted by MIT. The less severe impact scenario predicted fewer severed and heavily damaged 
columns when compared with the MIT and WAI studies. 

One reason for the higher damage predictions from WAI may be attributed to the lack of internal tower 
contents in the WAI model, such as workstations and other building contents. Another reason for the 
higher damage prediction in the WAI study could result from the aircraft model. In conducting a collapse 
analysis, WAI used engineering estimates to reduce the number of failed columns from that predicted by 
their FLEX model to stabilize the tower immediately after impact. Despite this adjustment, the WAI 
study still estimated significantly higher damage for WTC 1 than the MIT and NIST studies. For WTC 2 
their adjusted estimate fell in line with the MIT and NIST studies. 

11.10 SUPPORTING SIMPLIFIED ANALYSES 

Approximate analyses were carried out to provide guidance to the global finite element impact analyses. 
The specific analyses included: 

• The analysis of the overall aircraft impact forces: This analysis indicated that that the 
momentum transfer was more important than the strength of the impacting aircraft in 
determining the overall impact loads on the building. 

• The effects of the energy in the rotating engine components: This analysis indicated that the 
potential for the rotational kinetic energy of the engine to significantly change the magnitude 
of the tower impact damage is small. 

• The influence of the static preloads on the calculated impact damage: This analysis showed 
that the initialization of a static preload in the columns of the towers would not have a 
significant effect on the tower impact response, damage, or residual strength. 
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The analysis of the load characteristics required to damage core columns: This analysis 
determined the characteristics of a debris cloud (i.e. density and speed) necessary to fail a 
core column. Application of this analysis indicated that the dispersed fuel cloud alone would 
not be sufficient to fail core columns. 
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